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SUBJECT:  Transmittal  of  Technical  Report  D-78-34 


TO:  All  Report  Recipients 


1.  The  technical  report  transmitted  herewith  represents  the  results  of 
Work  Unit  1B08A  of  Task  IB,  Movements  of  Dredged  Material,  of  the  Corps 
of  Engineers'  Dredged  Material  Research  Program  (DMRP) . It  was  a part 
of  the  Environmental  Impacts  and  Criteria  Development  Project  (EICDP) , 
which  had  a general  objective  of  developing  techniques  for  determining 
the  spatial  and  temporal  distribution  of  dredged  material  discharged 
into  various  hydrologic  regimes.  The  study  reported  on  herein  was  part 
of  a series  of  research  contracts  developed  to  achieve  the  EICDP  general 
objective. 

2.  Regardless  of  the  location  or  character  of  a disposal  site,  an 
integral  part  of  the  problem  of  assessing  the  environmental  impact  of 
open-water  disposal  operations  is  an  ability  to  determine  the  movement 
of  dredged  material.  Dredged  material  at  an  aquatic  disposal  site  is 
subject  to  dispersion  by  the  tidal  stream,  estuarine  circulation,  waves, 
and  disturbances  of  the  hydraulic  flow  by  storms.  In  shallow  coastal 
waters,  sediment  may  be  transported  periodically  by  tidal  currents  or 
episodically  by  storm-generated  currents.  A major  objective  of  this 
investigation  was  to  evaluate  the  susceptibility  of  submerged  dredged 
material  mounds  to  these  disturbances. 

3.  This  report  describes  investigations  conducted  in  the  Gulf  of  Mexico 
at  the  Galveston,  Texas,  dredged  material  disposal  site  in  order  to 
determine  the  critical  erosion  velocity,  shear  stress,  and  modes  of 
sediment  transport  for  four  sediment  mixtures.  The  four  different 
sediments  eroded  similarly  as  evidenced  by  a similar  critical  bed  shear. 
Laboratory  results  were  compared  with  field  data,  indicating  greatest 
bedload  erosion  on  the  northern  portions  of  the  site  with  movement  of 
the  sediment  away  from  the  Galveston  Bay  entrance  channel. 

4.  The  information  and  data  published  in  this  report  are  contributions 
to  the  further  understanding  of  the  complex  nature  of  sediment  transport 
and  stability  of  submerged  dredged  material  deposits  and  should  help 
establish  a baseline  from  which  to  develop  meaningful  evaluations  for 
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the  selection  of  an  environmentally  compatible  disposal  alternative.  It 
s expected  that  the  interpretation  of  the  physical  interactions  will  be 
° significant  vaiue  to  those  persons  concerned  with  CE  dredged  material 
permit  programs. 


/JOHN  L.  CANNON 
Colonel,  Corps  of  Engineers 
Commander  and  Director 
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SUMMARY 


Flume  experiments  were  performed  using  four  sand,  silt, 
and  clay  mixtures  sampled  from  the  offshore  Galveston  dredged 
material  disposal  site  in  order  to  determine  the  critical  erosion 
velocity,  shear  stress,  and  modes  of  sediment  transport  for  each 
mixture.  Also,  an  analysis  of  the  hydrographic  regime  for  offshore 
Galveston  was  performed  based  on  meteoro logic  and  oceanographic  data 
collected  between  February  1975  and  June  1976.  Results  of  the  flume 
experiments  and  hydrographic  analysis  were  extrapolated  to  sediment 
transport  processes  believed  operative  in  the  offshore  disposal  site. 

The  results  of  the  flume  experiments  indicated  that  the 
four  sediment  mixtures  erode  similarly  as  evidenced  by  the  same  crit- 
ical bed  shear  (t^^.  ^ 1.0  dynes/cm^)  necessary  to  induce  massive 
bed  load  transport  and  rapid  suspended  sediment  increases.  Signifi- 
cant deposition  of  suspended  fine  silt  and  clay  occurred  at  low  bed 
shears  (x  < 0.2  dynes/cm"- ) . extrapolation  of  these  erosion  and 
deposition  bed  shears  to  recorded  offshore  bottom  current  speed 
measurements  over  month-long  periods  suggests  that  bed  load  erosion 
occurs  much  more  frequently  near  the  northern  margin  of  the  disposal 
site  than  at  the  southern  end.  In  addition , progressive  vector 
diagrams  Indicate  a net  down-coast  and/or  offshore-directed  bottom 
drift.  These  net  drift  directions  presumably  result  from  strong 
prevailing  onshore  winds  that  generate  onshore-directed  flow  near 
the  surface  and  down-coast  to  offshore  flow  at  the  bottom.  Maximum 
bottom  velocities  occur  during  these  wind-generated  bottom  flows. 


This  suggests  that  net  bedload  transport  of  disposal  mound  sediment 
is  oriented  away  from  the  Galveston  Bay  entrance  channel,  and  there- 
fore material  will  not  likely  return  as  shoaling  sediment. 
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PREFACE 


The  work  described  in  this  report  was  performed  under 
Contract  No.  DACW39-76-C-0I 1 5,  dated  30  June  IR7t>,  between  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg, 
Miss.,  and  the  Texas  A&M  University  (TAMED  Research  Foundation, 
College  Station,  Tex.  The  research  was  sponsored  by  the  Office, 
Chief  of  Engineers,  U.  S.  Army,  as  part  of  the  Dredged  Material 
Research  Program  (DMRP)  managed  by  the  Environmental  Laboratory 
(EL),  WES.  The  study  constituted  DMRP  Work  Unit  No.  IB08A. 

The  research  and  report  writing  were  conducted  by  and 
under  the  supervision  of  Anthony  J.  Moherek,  Department  of  Oceano- 
graphy, TAMU,  College  Station,  Tex.  Ms.  Jo  Ann  Treat  coordinated 
the  contract  for  the  TAMU  Research  Foundation.  The  hydrographic 
and  meteorologic  data  Incorporated  in  this  report  were  acquired 
from  work  conducted  under  Contract  No.  DACWe»4-75-C-0069. 

The  contract  was  managed  at  WES  by  Mr.  Barry  Holliday 
under  the  general  supervision  of  Dr.  Robert  M.  Engler,  Project 
Manager  for  the  Env i ronmenta I Impacts  and  Criteria  Development 
Project,  and  Dr.  John  Harrison,  Chief  of  EL. 


4 


Directors  of  WLS  during  the  study  and  preparation  of  this 
report  were  COL  G.  H.  Hilt,  CE,  and  COL  J.  L.  Cannon,  CE.  Technical 
Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 

UNITS  OF  MEASUREMENT 

U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (SI)  units  as  follows: 


Multiply 


*1 


To  Obtain 


feet 

miles  (U.  S.  statute) 
miles  (U.  S.  naut ica I ) 
i nch-ounces  ( force ) 
horsepower  (550  ft-lbf/s) 
gal  ions  (U.  S.  liquid) 
degrees  (angular) 


0.3048 
I .609344 
I .852 
0.0070615 
745.6999 
3.785412 
0.01745329 


metres 
kl lometres 
kl lometres 
newton-metres 
watts 

cubic  decimetres 
radians 


f5 


A rectanoular  circulating  tlumo  (Figure  I)  was  designed 


and  constructed  to  model  the  effects  of  unidirectional  turbulent  flow 


upon  bedload  and  suspended  sand,  silt,  and  clay  mixtures  sampled  from 


the  Galveston  dredged  material  disposal  site 


Figure  I.  Side  view  of  laboratory  flume 


The  critical  erosion  velocity,  in  terms  of  bed 


shear  and  current  speed,  necessary  to  erode  each 


sediment  mixture 


The  transport  conditions  during  which  this  eroded 


material  remained  in  suspension 


c.  The  amount  and  composition  of  suspended  matter 
that  resulted  from  erosion  under  controlled 
velocity  conditions. 

d . The  relationships  that  describe  the  recorded 
events  of  sediment  erosion  and  transport  in 
order  to  predict  the  relative  degree  of  erosion, 
transport,  and  deposition  occurring  in  the  off- 
shore dredged  material  disposal  site. 

In  addition,  this  study  included  an  evaluation  of  hydrographic  data 
acquired  within  and  Immediately  adjacent  to  the  offshore  disposal  site. 

These  data,  obtained  between  February  197b  and  June  l^7o,  constituted 
part  of  a related  study  (Estes  and  Scrudato  1977)  that  attempted  to 
determine  the  fate  of  dredged  material  discharged  into  the  Galveston 

site. 


Study  Area 

3.  The  location  of  the  offshore  disposal  site  is  presented 
In  Figure  2.  The  northern  boundary  is  situated  approximately  4.5  km 
southeast  of  the  southern  extremity  of  the  Galveston  Bay  entrance 
channel  south  jetty.  The  site  Is  near  rectangular,  extending  fo.7  km 
In  length  and  3.8  km  In  width,  and  Is  oriented  with  the  long  axis 
trending  In  a northwest-southeast  direction.  Also  shown  In  Figure  2 
are  the  locations  of  three  buoys  (B,  C,  and  D)  that  mark  specific 
dredged  material  discharge  sites  and  two  reference  areas  (Block  15  and 
Block  27)  that  served  as  control  sites  for  sediment  sampling  operations. 
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The  disposal  site  was  divided  Into  28  Individual  I ,3-km-square  blocks 
in  order  to  facilitate  sample  location  designations  (Figure  3). 

4.  Water  depths  within  the  disposal  site  ranged  from  10  m 
at  the  northern  extremity  to  a maximum  of  16  m at  the  southern  boun- 
dary. Prior  to  the  start  of  discharge  operations,  a bathymetric 
survey  (Figure  4)  was  performed  in  April  1975.  This  baseline  study 
Indicated  that  the  seafloor  of  the  disposal  site  is  very  smooth  with 
a dip  of  approximately  1.5  m per  mile*  towards  the  southeast.  Subse- 
quent to  this  survey,  dredged  material  discharge  operations  were 
carried  out  at  each  of  the  three  designated  discharge  sites  (Buoys  B, 
C,  and  D).  Iwo  bathymetric  surveys  have  been  completed  since  disposal 
operations  were  initiated  in  order  that  changes  in  the  disposal  mound 
topographies  over  time  due  to  sediment  transport  processes  could  be 
assessed  (Estes  and  Scrudato  1977).  The  bathymetric  configurations  of 
the  three  individual  disposal  mounds,  as  determined  by  the  June  1976 
survey,  are  illustrated  in  Figures  5,  6,  and  7. 


*A  table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (SI)  units  is  presented  on  page  15. 
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Figure  6,  Bathymetric  configuration  of  sediment  mound  at  Buoy  C,  Contour  interva 
(from  Estes  and  Scrudato  1377) 


Bathymetric  con* igurat ion  of  sediment  mound  at  Buoy  D,  Contour  interval 
(from  Estes  and  Scrudato  1977) 


PART  II;  LITERATURE  REVIEW 


General 

5.  The  physical  parameters  governing  sediment  transport 
processes  of  marine  cohesive  material  are  at  best  qualitatively 
understood.  This  is  due  to  the  fact  that  the  focus  of  recent  hydra- 
ulic research  has  centered  on  the  erosion,  transport,  and  deposition 
of  primarily  loose,  noncohesive  particles  such  as  rounded  quartz 
grains.  To  the  author's  knowledge,  no  studies  have  been  reported 
dealing  directly  with  sediment  mixtures  containing  near  equivalent 
amounts  of  sand,  silt,  and  clay. 

6.  The  first  significant  study  establishing  erosion-deposition 
criteria  for  sedimentary  particles  of  uniform  size  was  presented  by 
Hjulstrom  (1935).  Figure  8 illustrates  the  well  known  "Hjulstorm 
diagram,"  The  flow  velocities  given  represent  the  average  current  speed 
that  is  presumably  40  percent  greater  than  the  bottom  velocity  tor  a 
flow  depth  exceeding  I m.  The  diagram  shows  that  erosion  of  loose, 
noncohesive  fine  sand  occurs  at  significantly  lower  current  velocities 
than  for  cohesive  silts  and  clays.  Alternately,  tine  particles  are  apt 
to  remain  in  suspension  at  considerably  lower  flow  speeds  than 

coarse  particles.  Postma  (1967)  presented  a refinement  of  Hjulstorm's 
diagram  using  data  reported  by  Sundborg  (1956).  He  related  sediment 
transport  criteria  to  varying  degrees  of  cohesive  sediment  consolida- 
tion. The  results  depicted  in  Figure  9 are  not  likely  to  be  represen- 
tative of  erosion-deposition  criteria  for  discharged  dredged  material 
having  highly  variable  consolidation  characteristics. 
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Bed  Shear  Stress 


7.  Workers  realized  early  that  description  of  the  competency, 
or  erosion  capability,  of  currents  acting  on  a given  bedload  depended 
upon  an  accurate  measure  of  the  turbulent  flow  intensity.  For  this 
reason,  hydraulic  engineers  defined  the  term  bed  shear  stress  (t  )*  as 
representing  the  frictional  force  acting  upon  the  sediment-water  inter- 
face. The  critical  shear  stress  (t  ) was  then  defined  as  that  force 

ocr 

necessary  to  cause  incipient  bedload  erosion.  Various  empirical  tech- 
niques have  been  proposed  for  measuring  the  bed  shear  within  a given 
flow.  Bed  shear  for  open  channel  flow  has  been  commonly  expressed  as: 


where  (y)  is  defined  as  the  unit  weight  of  the  fluid;  (R^)  represents 
the  hydraulic  radius  or  wetted  perimeter  of  the  channel;  and  (S)  Is 
the  slope  of  the  energy  grade  line.  This  expression  is  representative 
of  the  average  value  of  the  shear  stress  per  unit  wetted  area.  However, 
the  product  (yR^S)  is  not  applicable  to  the  marine  environment  since 
the  slope  of  the  energy  gradient  cannot,  at  this  time,  be  directly 
measured.  Thus,  most  Investigators  resorted  to  usage  of  Prandtl-von 
Karman's  universal  velocity  distribution  law  for  fully  developed 
turbulent  flow: 


U 


2.5 


I n 


(2) 


•For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  F), 
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where:  U 


the  flow  speed  at  elevation  (Y)  above  the  bed 
the  roughness  element 


I 


I 


Zo  « 

u*  = -^t0/p  = the  friction  element 

p = water  density 

8.  Simple  regression  methods  can  be  employed  to  estimate 

(u»)  and  (Z.)  in  order  to  calculate  (t  ).  The  main  restriction  of 
o o 

this  relationship  is  that  current  velocity  measurements  are  required 
from  at  least  two  points  immediately  above  the  bed.  Teteki  (1972) 
related  Von  Karman's  logarithmic  law  to  the  "velocity  defect  law"  and 
"law  of  the  wall"  in  turbulent  boundary  layers  of  incompressible  steady 
flow. 

The  Problem 


9.  A universal  relationship  equating  erosion  and  transport 
of  cohesive  sediment  has  yet  to  be  formulated  due  to  the  complex 
adhesive  and  cohesive  forces  acting  between  sediment  particles.  Graf 
(1971)  outlined  the  basic  problem  by  defining  the  critical  shear  stress 
necessary  to  erode  sediment  as  a function  of  sediment  particle  diameter 
(d)  and  a coefficient  of  cohesion  (CQ),  such  that: 


Tocr  - f(d'  co> 


(3) 


Previous  Flume  Experiments 

10.  Flume  experiments  performed  by  Smerdon  and  Beaslay  (1969) 
using  II  different  Missouri  soils  demonstrated  that  the  critical  shear 
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stress  correlated  best  with  the  plasticity  Index  (lp)  of  the  soil.  This 
study  used  Equation  I to  calculate  the  bed  shear,  or  unit  tractive  force, 
acting  on  the  soil  for  a given  flow  rate,  defining  incipient  bed  load 
motion  as  corresponding  to  "general  movement  of  the  bed  material."  Sim- 
ilar experiments  by  Lyle  and  Smerdon  (1965),  testing  seven  different 
Texas  soils,  determined  that  (lp)  and  the  void  ratio  (e)  correlated  well 
with  (t  ).  Determination  of  the  critical  shear  stress  was  made  by 
plots  of  the  soil's  erosion  rates  against  the  measured  bed  shear.  The 
point  at  which  the  erosion  rate  increased  sharply  indicated  initial  bed 
failure.  For  both  studies,  results  were  based  upon  erosional  current 
flow  over  a smooth  plane  bed.  The  applicability  of  their  results  to 
dredged  material  discharged  into  an  open  marine  environment  appears 
highly  tenuous  because  of  the  hydrodynamics liy  rough  surface  produced 
by  material  mounding  on  the  sea  floor. 

II.  Einstein  and  Krone  (1962)  correlated  modes  of  sediment 
transport  with  suspended  sediment  concentration  using  San  Francisco 
Bay  mud  containing  nearly  equal  amounts  of  silt  and  clay.  A recircu- 
lating flume,  1.0  m wide  and  30.48  m long,  was  filled  to  a 30.5-cm 
water  depth  for  monitoring  of  the  suspended  matter  concentrations  at 
uniform  flow  speeds.  Depositional  concentration  changes  indicated 
that  below  300  ppm  the  suspended  sediment  concentration  decreased 
logarithmically  with  time.  Plots  of  deposition  rates  versus  bed  shear 
showed  that  at  shears  exceeding  0.60  dyne  / cm 2 (kPa)  deposition  of 
suspended  particles  did  not  occur.  For  concentrations  greater  than 
300  ppm,  flocculation  of  the  suspended  clays  became  a dominant  factor 
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in  increasing  deposition  rates.  This  was  attributed  to  "an  increase 
in  frequency  of  interparticle  contacts  forming  large  floe  sizes  over 
a short  time  period."  Scour  experiments,  performed  by  Krone  (1962) 
with  the  same  bed  load  material,  revealed  a log- log  relationship  between 
the  suspended  sediment  concentration  and  time,  suggesting  that  exchange 
of  suspended  and  deposited  sediment  occurred.  No  correlation  was  found 
between  the  critical  shear  stress  and  the  shear  strength  of  the  bed  as 
determined  by  in  situ  penetrometer  measurements. 

12.  Partheniades  and  Mehta  (1971)  presented  data  on  rates  of 
deposition  of  fine  cohesive  sediments  in  turbulent  flows.  Using  kaolinite 
clay  suspensions  in  a rotating  annular  channel  and  ring  apparatus,  they 
demonstrated  that  the  percentage  of  the  total  sediment  within  a given  flow 

depends  only  on  the  bed  shear  stress  (t  ).  Furthermore,  the  percentage 
(C1)  of  the  depositable  sediment  deposited  at  time  (t)  was  found  to  vary 
with  time  according  to  the  law. 


C'  * a log  t + 6 (4) 

where  a is  a constant  independent  of  the  flow  conditions  and  sediment 
concentration,  and  0 is  a constant  dependent  on  shear  stress  only. 

Offshore  Sediment  Transport 


13.  Understanding  of  the  extent  of  sediment  transport  processes 

in  the  open  marine  environment  requires  knowledge  of  the  fluid  motions 

< 

involved  both  in  the  entire  water  column  and  at  the  sediment-water 
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interface.  Weggel  (1972)  outlined  the  fundamental  fluid  motions  in 
the  water  column,  as  well  as  their  respective  generating  forces.  These 
forces  Include:  meterologlcal , such  as  wind-driven  surface  stresses 
generating  waves  and  storm  surges;  astronomical  forces  producing  tides 
and  associated  tidal  currents;  and  baroclinic  forces  induced  by  fluid 
density  gradients  that  are.  In  turn,  due  to  temperature  and  salinity 
differences.  Figure  10  illustrates  schematically  the  relationship 
between  these  forces  and  flow  regions  within  the  water  column  (0  is 
mean  flow  speed). 

14.  Telekl  (1972)  described  the  complexities  in  evaluating 
near-bottom  oscillatory  motion  generated  by  surface  waves.  His  study 
demonstrated  that  sediment  motion  is  dependent  upon  the  characteristic 
frequency  and  amplitude  of  oscillation  as  determined  experimentally  by 
shear  stress  measurements  shoreward  outside  the  zone  of  breaking  waves. 
Komar  and  Miller  (1973)  derived  an  empirical  relationship  for  determin- 
ing the  oscillatory  bottom  velocities  needed  to  reach  the  threshold  of 
sediment  movement  for  grain  diameters  less  than  0.05  cm  (medium  sand). 
Phillips  (1969)  pointed  out  that  surface  waves  induce  a mean  forward 
velocity  near  the  bottom  in  addition  to  strictly  oscillatory  motion. 
This  phenomena,  first  observed  by  Caligny  (1878),  has  been  termed 
boundary  layer  streaming.  During  the  literature  survey,  no  studies 
were  found  that  quantitatively  assessed  the  effects  of  near-bottom 
oscillatory  motions  on  heterogeneous  sediment. 

15.  The  theory  of  surface  wind  stresses  generating  nearshore 
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STORM  SURGES 


10.  Relationship  between  generating  forces  and  flow  regions  within  the  offshore  water  column 


vertical  circulation  patterns  has  been  discussed  by  Sverdrup  et  al., 
(1942).  Murray  (1972)  presented  results  of  current  and  water  density 
observations  recovered  at  an  anchor  station  in  15  m of  water  east  of 
the  Mississippi  River  delta.  Measurement  of  the  currents  plus  salinity- 
induced  density  gradients  for  a 5-day  period  in  March  1970  showed  that 
southeasterly  storm  winds  exceeding  15  m/sec  produced  a circulation 
pattern  causing  onshore  flow  near  the  surface  and  offshore  flow  at 
the  bottom.  Currents  associated  with  the  storm  persisted  for  at  least 
3 days  after  the  storm  passage,  generating  stable  offshore  bottom 
currents.  Longard  and  Banks  (1952)  also  reported  that  wind-driven, 
shoreward-directed  surface  waters  were  coupled  to  of f shore-moving  bottom 
currents  at  a station  in  91  m of  water,  southeast  of  Halifax,  Nova 
Scotia. 


Flow  Near  the  Bed 

16.  Criteria  for  establishing  whether  a given  flow  is  laminar 
or  turbulent  in  nature  is  based  upon  the  calculated  value  of  its 
Reynolds  number  (Re)  such  that: 

Re  = HtP  (5) 

e u 

where  (L)  is  the  reference  length  and  (p)  represents  dynamic  viscosity 
(after  Reynolds  1883).  The  critical  value  of  the  Reynolds  number  for 
the  transition  from  laminar  to  turbulent  flow  depends  upon  the  repre- 
sentative length  scale,  velocity,  and  geometry  of  the  flow  syste  i. 

Based  on  numerous  experiments  with  flow  through  pipes,  Reynolds  deter- 
mined that  flows  with  (R0)  <2000  are  laminar,  and  those  characterized 
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by  (Re)  >2000  are  turbulent.  For  oceanic  currents,  it  is  apparent 
that  the  representative  length  scale  will  be  large  resulting  in  Reynolds 
numbers  greater  than  a few  thousand  for  current  speeds  exceeding  a few 
centimeters  per  second.  Thus,  such  currents,  on  a macroscopic  scale, 
can  be  considered  as  fully  turbulent. 

17.  Prandtl  and  Tietzens  (1934)  demonstrated  that  the  turbu- 
lent layer  immediately  above  the  bed  consisted  of  a flow  region  where 
the  horizontal  component  of  the  fluid  velocity  approaches  zero.  This 
zone  has  commonly  been  observed  as  being  laminar  (or  streamline)  and  is 
referred  to  as  the  laminar  (or  viscous)  sublayer.  Figure  II  illustrates 
the  velocity  distribution  in  the  turbulent  boundary  layer  (modified  from 
Rouse  1961). 

18.  Nikuradse  (1933)  demonstrated  that  hydraulically  rough 
or  turbulent  flow  occurs  when  the  scale  length  of  the  bottom  roughness 
(e.g.,  median  grain-size  diameter)  exceeds  the  laminar  sublayer  thick- 
ness. Under  such  conditions,  irregularities  on  the  bed  cause  local 
variations  in  the  bed  shear  stress  and  consequently  reduce  the  average 
bottom  shear  stress  needed  to  initiate  grain  movement.  For  inhomo- 
geneous discharged  dredged  material  on  the  sea  floor,  it  seems  appro- 
priate to  classify  flow  near  the  bed  as  hydrodynami ca I I y rough  resulting 
in  locally  wide  variations  in  the  bed  shear  stress.  The  implication  of 
this  is  that  determination  of  the  critical  shear  stress  needed  to  cause 
incipient  bed  load  erosion  becomes  difficult.  Additionally,  bottom 
sediment  is  likely  to  be  entrained  within  a given  flow  above  a rough 
bed  at  lesser  current  speeds  than  for  a smooth  bed. 
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PART  III:  TECHNIQUES 


Field  Techniques 

Offshore  sediment  sampling 

19.  Offshore  sediment  sampling  operations  were  carried 
out  onboard  the  R/V  TEXAS  STAR. 

20.  A 0.3048-cu-m  box  core  was  used  to  recover  approxi- 
mately 100  kg  of  sediment  from  four  localities  within  the  field 
study  area.  These  sampling  stations  were  designated  as  Block  15, 

Block  27,  Buoy  C,  and  Buoy  D.  Sediment  samples  from  Block  15  and 
Buoy  C were  obtained  during  November  1975;  Block  27  and  Buoy  D sites 
were  sampled  In  May  1976.  Blocks  15  and  27  were  selected  as  reference 
sites  since  they  were  believed  to  be  uncontaminated  by  dredged  material 
disposal  operations  being  carried  out  at  Buoys  B,  C,  and  D.  Figure  3 
depicts  the  locations  of  the  sampling  stations.  An  inoperative  Del 
Norte  navigation  system  on  board  the  R/V  TEXAS  STAR  prevented  accurate 
location  of  the  reference  blocks.  Thus,  the  Block  15  station  was 
located  by  dead  reckoning  from  Buoy  B.  Similarly,  the  position  of 

the  Block  27  station  was  ascertained  by  dead  reckoning  from  Buoy  D. 
Buoy  C and  D sampling  stations  were  located  approximately  30  m south- 
west of  the  respective  buoys.  Such  positioning  permitted  sampling  of 
previously  deposited  dredged  material.  All  sediment  samples  were 
stored  in  sealed  55-gal  drums  until  experiments  could  be  performed. 

Collection  of  current  data 

21.  Two  types  of  current  meters  were  employed  in  the  field 
operations:  a hand-held  Bendix  Q-15  current  meter  and  two  Bralncon 
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Savon i us  continuous-recording  meters.  Malfunctioning  of  the  current 
meters  did  occur  occasionally  during  scheduled  periods  of  observations. 
Table  I contains  exact  time,  location,  and  type  of  current  meter  data 
recovered . 


22.  Vertical  profiles  were  taken  with  the  0-15  meter  at  Buoy 
B,  C,  and  D stations.  The  profiles  were  obtained  by  lowering  the  Q-15 
meter  from  the  surface  In  1.5-m  intervals  to  a position  1.5  m above  the 
bottom.  Then  It  was  lowered  in  0.3-m  increments  to  the  seabed.  In 
this  manner,  point  measurements  of  the  ambient  current  speed  and  direc- 
tion were  recorded  for  the  entire  water  column  with  emphasis  on  the 
lowermost  1.5  m where  profiles  of  the  bottom  currents  were  desired, 

23.  Bottom  current  measurements,  obtained  with  Savonius 
current  meters,  consisted  of  current  speed  and  direction  recorded  at 
1 5-mi n intervals,  1.0  m above  the  bottom.  Bottom  current  speed  and 
azimuth  were  recorded  for  month-long  periods  at  Buoy  B and  0 stations. 
Figure  12  schematically  illustrates  the  bottom  current  meter  arrays 
positioned  approximately  50  m southwest  of  Buoys  B and  D.  A submer- 
sible buoy  above  and  approximately  175  kg  of  anchor  weight  below  held 
This  instrument  stationary  1.0  m off  the  bottom. 

Geotechnical  Properties  of  Bed  load 

24.  The  geotechnical  properties  measured  on  the  field  samples 
obtained  from  the  Galveston  site  include  (I)  grain-size  distribution 
(dispersed  and  nondlspersed  analysis),  (2)  water  content,  (3)  bulk 
density,  (4)  void  ratio,  (5)  specific  gravity,  (6)  remolded  shear 
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Table  I 


Type,  Location.  Date,  and  Time 
of  Current  Meter  Data  Recovered 


TVPe 

Velocity 
Prof  1 les 


Bottom 
Current 
Read  I ngs 


Location 

Date 

Time  (CST) 

Buoy  B 

8 April  1976 

1000,  1120  hr 

9 April  1976 

1055  hr 

1 May  1976 

1150  hr 

Buoy  C 

8 Aprl 1 1976 

1430  hr 

1-2  May  1976 

1400,  1600,  1800,  2000,  2200, 
2400,  0200,  0400,  0600,  0800, 
1000.  1200.  1400  hr 

Buoy  D 

8 April  1976 

1245,  1330  hr 

3 May  1976 

1200  hr 

Buoy  B 

17-24  March  1976 

1100-1315  hr 

18  Aprl 1 - 1 May 
1976 

1200-0915  hr 

1 May  - 1 1 May 

1976 

1345  to  1930  hr 

Buoy  D 

9 Oct.  - 14  Nov. 
1975 

1600  - 1020  hr  (from  Hall 

1976) 

17  March  - 8 

April  1976 

1230  - 1030  hr 

8 Apr  I I - 3 May  1415  - 0900  hr 
1976 


3 May  - II  May 
1976 
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1345  - 0715  hr 


Figure  12.  Bottom  current  meter  array  positioned  adjacent  to  both 
Buoy  B and  Buoy  D 

40 


L 


strength,  (7)  organic  content,  and  (8)  Atterberg  limits  (liquid  limit, 
plastic  limit,  and  plasticity  index).  Each  analysis  was  performed  on 
homogenized  sediment  samples  obtained  just  prior  (within  minutes)  to 
their  introduction  into  the  flume.  Homogenization  entailed  thorough 
mixing  of  the  sample  with  a stirring  rod  until  ail  clay  lumps  were 
dispersed.  The  following  sort  I on  describes  the  individual  techniques 
that  were  employed  for  each  geotechnical  property  analysis. 

Textural  analysis 

25.  The  dispersed  grain-size  analysis  performed  on  the 
homogenized  field  samples  followed  the  procedure  given  by  Folk  (1974). 

A solution  of  Calgon  and  distilled  water  was  added  to  approximately 

10  g of  dredged  material.  This  mixture  was  then  shaken  to  disassociate 
the  clay  particles  before  being  sieved  through  a U.S.  Standard  No.  230 
sieve  (0.0625-mm  openings)  into  a I cylinder.  The  filtrate  was 
pipetted  according  to  Folk's  method.  The  material  retained  on  the  230 
mesh  sieve  was  then  dried  and  sieved  through  sieve  sizes  of  U.S. 

Standard  No.  10  (2.0  mm),  18  (1.0  mm),  35  (0.50  mm),  60  (0.25  mm),  IZ0 
(0.125  mm),  and  230.  Finally,  the  sediment  was  weighed  using  a precision 
balance  (accuracy  0.0005  g)  in  order  to  calculate  the  proportion  of 
material  per  given  grain-size  ($)  size  range.  This  proportion  was 
recorded  as  percent  dry  weight. 

26.  The  nondispersed  grain-size  analysis  followed  the  pro- 
cedure cited  above  except  that  distilled  water  was  substituted  for 
Calgon  dispersing  solution  when  called  for.  This  physical  condition 


] 
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was  desired  in  order  to  assess  the  grain-size  distribution  of  the  field 
samples  in  a manner  that  more  closely  approximated  their  In  situ  state. 

27.  Lastly,  cumulative  curves  were  plotted  from  the  cumula- 
tive weight  percent  data  for  each  phi  size.  Since  the  No.  10. 2-^ 
size  range  was  the  highest  phi  size  measurement  (smallest  grain  size), 
that  part  of  the  cumulative  curve  was  extrapolated  to  higher  phi  sizes 
(>10.2  ) by  continuing  the  curve  upward  at  approximately  the  equivalent 
rate  of  the  slope  chanae  determined  by  the  9-  and  10.2-0-size  data 
points.  Each  cumulative  curve  supplied  the  necessary  information  to 
calculate  statistical  parameters  that  described  the  grain-size  distri- 
bution. The  parameters  calculated  were  graphic  mean,  inclusive  graphic 
standard  deviation,  and  inclusive  graphic  skewness.  These  parameters, 
from  Folk  (1974),  are  as  follows: 

= Graphic  Mean 

= 016  + $50  + 084  (6) 

3 


0j  = Inclusive  Graphic  Standard  Deviation 

= 084  - 016  + p95  - p5 

4 §76 


( 7 ) 


SKj  = Inclusive  Graphic  Skewness 


= Q 1 6 t Q84  - 2 Q50  + Q5 


2(084  - 016) 


t p95  - 2 060 

(095  - 08) 


(8) 


Shear  strength 

28.  The  vane  shear  test  consisted  of  Immersing  a vane  shear 
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into  a 50-g  homogenous  sample  and  measuring  the  torque  required  to 
shear  the  sediment  around  the  vane's  cylinder  of  rotation.  The  low 
shear  strength  of  the  field  samples  necessitated  use  of  a visco-elastic 
vane  shear  for  sample  measurements.  This  instrument,  developed  by 
Stevenson  (1973)  for  determination  of  the  visco-elastic  modulus  of  sub- 
marine sediments,  consists  of  a four-btaded  vane-shear  coupled  to  a 
variable-speed  motor  via  a step-function  reducer. 

29.  Torque  was  measured  by  a strain-gaged  torque  transducer 
and  was  recorded  on  an  oscillograph  as  a plot  of  torque  versus  rota- 
tion angle.  Torque  sensitivity  of  this  instrument  was  measured  to  be 

_5 

1.38  X 10  kq-m  (0.001  in. -lb).  Point  of  maximum  torque,  termed 
failure,  was  read  from  the  torque  versus  rotation  angle  plot  where 
torque  values  became  constant  with  increasing  rotation  angle.  Figure 
13  represents  an  example  of  this  plot.  Conversion  of  the  maximum 
torque  to  standard  units  of  shear  strength,  kilopascals  (kPa),  was 
accomplished  by  simple  multiplication  of  torque  units  (inch-ounces)  by 
the  factor  1.3434  X 10“^. 

Void  ratio,  water  content, 

bulk  density,  and  specific  gravity 

30.  Void  ratio,  water  content,  bulk  density,  and  specific 
gravity  were  determined  from  a single  sample.  Determination  of  these 
geotechnical  properties  generally  followed  the  methodology  presented 
by  Cernock  (1967).  Correction  for  salt  content  in  the  interstitial 
water  of  the  sediment  was  excluded  because  (I)  measurements  of  the 
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Figure  13.  Example  of  torque  versus  rotation  angle  plot 


I 


Interstitial  water  salinity  were  lacking,  and  (2)  the  weight  contri- 
bution of  the  salts  to  the  sediment  mass  was  negligible,  since  the 
field  area  bottom  water  salinity  was  less  than  28  °/oo  (from  Cool  1976). 
Prior  to  Introduction  Into  the  flume,  the  test  material  was  sampled 
twice  In  order  that  two  values  of  each  physical  property  could  be  ascer- 
tained. These  two  values  were  then  averaged. 


31.  The  laboratory  procedure  for  this  analysis  required  use 
of  the  following  apparatus:  (I)  a Beckman  Air  Comparison  Pycnometer 
utilizing  helium  gas  for  volumetric  determinations  of  sediment  samples 
(accuracy  - 0.02  cm^),  (2)  a Sartorius  Electra  balance  for  weight 
measurements  (accuracy  - 0.0001  g),  and  (3)  a Thelco  oven  for  sample 
dessicatlon  at  I05°C. 


32.  The  analytical  procedure  consisted  of  first  placing 
approximately  40  g of  homogenized  sediment  Into  a tared  beaker  for 
Immediate  wet  volume  and  wet  weight  measurement.  Next,  the  sample  was 
oven-dried  at  I05°C  for  24  hr  and  then  allowed  to  equilibrate  for  I to 
2 hr  In  a dessicator.  The  dry  weight  of  the  sample  was  obtained  prior 
to  mechanical  crushing  via  mortar  and  pestle  for  volumetric  determina- 
tion of  the  solid  mass.  An  additional  24-hr  oven-drying  period  after 
crushing  removed  the  interstitial  moisture. 


33.  Final  recording  of  the  crushed  dry  volume  and  dry  weight 
permitted  computation  of  the  four  previously  cited  physical  properties. 
The  equations  utilized  for  these  calculations  are  as  follows  (modified 
after  Cernock  1967): 

a.  Void  ratio  (e)  is  defined  as  the  ratio  of  void  volume 


(Vv)  to  the  volume  of  the  solid  particles  (Vg)  within  a given  sediment 
mass,  such  that 

Vv 

e = VT  (9) 

where  (Vv)  is  defined  as  the  total  volume  of  the  sediment  (VT)  minus  the 
volume  of  the  particles  (Vs). 

il*.  Water  content  (W)  Is  defined  as  the  ratio  in  percent  of 
the  weight  of  seawater  (Wsw)  to  the  weight  of  the  oven-dried  sediment 
(Ws)  such  that 

W, 


sw 


Wc 


100 


(10) 


£jl  Bu I k dens  1 ty  (Y^)  is  defined  as  the  ratio  of  the  total 
weight  of  the  sediment  mass  (WT)  to  the  total  volume  of  the  sediment 
mass  (Vy)  such  that 

W, 


T 

Yh  - 7— 


(II) 


1*.  Specific  gravity  (G)  is  defined  as  the  ratio  of  weight  in 
air  of  a given  solid  volume  (Ws)  to  the  weight  in  air  of  an  eguivalent 
volume  of  distilled  water  at  4°C  (VSYW)  such  that 

Wc 


G = v Y 


s 

sYw 


(12? 


where  (Yw)  is  defined  as  the  specific  gravity  of  distilled  water  at 
4°C  (1.000). 

Atterberg  limits 

34.  The  procedure  for  the  determination  of  the  Atterberg 
limits  (i.e.,  liquid  limit,  plastic  limit,  and  plasticity  Index)  followed 
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that  given  by  Lambe  (1951).  Due  to  the  high  sand  content  10  percent) 
in  each  material  type,  all  samples  required  wet-sieving  using  distilled 
water  through  a U.S.  Standard  No.  40  (0.42-mm)  sieve  In  order  to  produce 
a test  sample  that  was  within  the  upper  limit  of  grain  sizes  ( 1.25$) 
used  for  Atterberg  limit  analysis. 

35.  Liquid  limit.  Liquid  limit  (W[_)  is  defined  as  the  mois- 
ture content,  percent  dry  weight,  of  the  oven-dried  sediment  at  which 
the  material  will  flow  to  close  a 2-mm-wide  groove  when  jarred.  Each 
material’s  liquid  limit  was  determined  by  first  subdividing  the  sieved 
fraction  into  four  equal  aliquots.  Next,  one  aliquot  at  a time  was 
homogenized  into  a uniform  paste  and  placed  into  the  cup  of  a liquid 
limit  device.  A special  grooving  tool  of  standard  dimensions  was  drawn 
through  the  sample  resulting  in  a groove  2-mm  wide.  The  cup  was  cali- 
brated 1.0  cm  above  a flat  surface  upon  which  it  was  mechanically  jarred 
when  operative.  The  number  of  blows  required  to  close  the  groove 
approximately  1.25  cm  along  its  length  was  recorded.  A 10-g  sample  was 
immediately  extracted  from  the  closed  groove.  This  sample  was  placed 
into  a tared  container  and  weighed  with  a Sartorius  Electra  balance  for 
wet  weight  determination.  This  process  was  repeated  for  each  aliquot. 
Finally  all  four  tested  aliquots,  per  given  sediment  type,  were  oven- 
dried  for  24  hours  and  then  weighed.  Knowing  the  dry  and  wet  weight 
of  each  aliquot  permitted  calculation  of  its  moisture  content.  Then, 
a plot  of  the  moisture  content  versus  thQ  log  of  the  number  of  blows 
required  to  close  each  groove  was  drawn.  Visual  interpolation  of  a 
regression  line  through  the  four  plotted  points  of  measurement  allowed 
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computation  of  the  liquid  limit.  This  value  corresponded  to  the 


moisture  content  at  log  25  blows  (after  Lambe  1951). 

36.  PI ast ic  limit.  The  plastic  limit  (Wp)  is  defined  as  the 
lowest  moisture  content  of  the  oven-dried  sediment  at  which  the  material 
can  be  rolled  into  threads  0.318  cm  in  diameter  without  breaking  into 
pieces.  After  homogenizing  approximately  15  o of  the  material,  the 
sample  was  hand-rolled  on  a smooth  glass  plate  continuously  until  a 
thread  0.318  cm  in  diameter  formed  and  crumbling  appeared.  This  thread 
was  placed  in  a tared  container  for  water  content  determination  according 
to  the  method  described  in  the  previous  section.  This  moisture  content 
was  then  recorded  as  the  plastic  limit. 

37.  Plasticity  index.  The  plasticity  index  (lp)  is  defined  as 
the  difference  between  the  liquid  and  plastic  limit.  Hence,  determination 
of  the  latter  properties  permitted  simple  computation  of  (Ip). 

38.  Organic  content.  Six  dried  samples  analyzed  for  combus- 
tible organic  content  were  placed  In  a Themco  infrared  oven  and  heated 
for  2 hr  at  600°C.  Measurement  of  the  weight  difference  before  and 
after  combustion  permitted  computation  of  the  combustible  organic  con- 
tent as  percent  dry  weight.  This  process  was  considered  an  approxi- 
mation since  inorganic  particulate  matter  was  partially  oxidized  in 
the  process.* 


•Personal  communication,  David  B.  Mathis,  U.  S.  Army  Engineer  Waterways, 
Experiment  Station,  3 May,  1976. 
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Geotechnical  Properties  of  Wash  load 


39.  One-half-liter  water  samples  were  collected  at  specific 
time  Intervals  (see  Appendix  A)  from  the  flume  during  experimental  runs 
in  order  to  ascertain  total  suspended  matter  concentration,  suspended 
sediment  grain-size  distribution,  and  washload  composition.  Emphasis 
was  placed  upon  measuring  the  total  suspended  matter  concentration  of 
each  sample  obtained  since  changes  in  the  wash  load  concentration  were 
expected  to  be  a function  of  bed  load  erodibility  or  stream  competency. 
The  following  sections  describe  the  laboratory  techniques  employed  to 
quantitatively  assess  the  physical  characteristics  of  the  washload. 

Total  suspended  matter  concentration 

40.  The  total  suspended  matter  concentration  (TSM)  of  the 
washload  samples  was  determined  by  the  gravimetric  procedure  outlined 
by  Cool  (1976).  Individual  47-mm-diam  Nucleopore  GE-40  filters  (pore 
size  = 0.40  pm)  were  each  preweiqhed  twice  using  an  Ainsworth  Type  24  N 
balance  (accuracy  0.0001  g).  Expected  static  electricity  was  minimized 
during  the  weighing  process  by  passing  each  filter  over  uranyl  acetate 
crystals  and  then  weighing  each  one  beneath  an  a-emitting  ionizing 
source  (Polonium  210).  Three  control  filters  were  each  weighed  five 
times  during  the  weighing  of  each  set  of  experimental  filters.  Thus, 
a correction  factor  was  attained  to  compensate  for  changes  due  to 
humidity,  temperature,  etc.,  between  the  time  the  filters  were 
preweighed,  used,  and  dried.  Control  filter  weight  variations  were, 
found  to  be  approximately  8 to  10  pg,  which  was  considered  negligible 
compared  to  the  concentration  levels  measured  during  the  experimental 
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runs 


41.  Apparatus  utilized  for  the  suspended  matter  filtration 
was  composed  of  five  mil  I i pore  in-line  disc  filter  holders,  0.5-8. 
plastic  sample  containers,  Tygon  tubing,  one  Gast  760-mm  mercury 
vacuum  pump,  and  five  I .0-i-vol ume  vacuum  flasks.  Figure  14  illus- 
trates schematically  the  filter  system  used. 

42.  Double-distilled  deionized  water  was  used  to  rinse  and 
thus  dissolve  retained  sea  salts  on  the  Nucleopore  filters.  The 
filtrate  was  tested  with  0.1  silver  nitrate  for  evidence  of  salt 
contamination.  Three  to  four  thorough  washings,  approximately  25  ml 
each,  removed  all  visual  traces  of  salt  precipitation  If  the  filters 
were  allowed  to  drain  dry  before  the  next  wash.  All  filters  were 
stored  in  a desslcator  for  at  least  5 days  prior  to  weighing  to  ensure 
thorough  drying.  The  total  suspended  sediment  concentration  in  mg  per  e 
was  calculated  by  dividing  the  weight  difference  of  the  dried  filters 
by  the  volume  of  water  filtered. 

Grain-size  distribution 

43.  A Model  TA  Coulter  Counter  was  employed  for  determining 
the  suspended  sediment  grain-size  distribution  of  selected  water  samples 
taken  during  experimental  runs  with  the  flume.  The  laboratory  procedure 
given  by  Shldeler  (1975)  with  certain  modifications  was  followed. 

Sheldon  and  Parsons  (1966),  Shldeler  (1975),  and  others  have  established 
through  repetitive  experimentation  that  the  size  distribution  of  a 
particular  sample  can  be  determined  using  a Coultor  Counter  with  a 
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precision  equivalent  to  or  better  than  the  standard  time-consuming 
pipette  analysis. 

44.  During  each  analysis,  two  aperture  tubes  with  200-um 
and  30-um  apertures,  respectively,  were  utilized  to  provide  a total 
analytical  size  range  of  64  pm  to  0.630  pm  (4.0  to  1 0 . 6<J> ) . Routine 
mud  analysis  performed  by  Shideler  (1975)  found  that  the  combination 
of  the  two  tubes  was  optimal  in  deciphering  suspended  matter  grain- 
size  distributions. 

45.  Since  each  tube  was  only  capable  of  measuring  particles 
in  a size  range  from  2 to  40  percent  of  the  tube  aperture  diameter, 
each  sample  required  appropriate  sieving  before  being  analyzed.  Thus, 
each  sample  was  wet-sieved  through  a U.S.  Standard  No.  120  (125-pm) 
sieve  prior  to  performing  the  200-pm  aperture  tube  analysis.  This 
sample  was  then  wet-sieved  through  a 20-pm  micro-mesh  sieve  for  the 
30-pm  aperture  tube  analysis.  After  recording  the  size  distribution, 
percent  total  volume,  measured  by  each  aperture  tube,  the  "two-tuba' 
results  were  arithmetically  combined  and  normalized  to  100  percent 
volume  distribution  according  to  the  overlap  procedure  given  by  Coulter 
Electronics  (1975).  Figure  15  Illustrates  a typical  worksheet  with 
recorded  results  of  one  analysis. 

46.  After  data  normalization,  the  resultant  cumulative  volume 
percent  for  each  size  range  was  plotted  on  probability  scale  graph  paper 
to  produce  a cumulative  curve.  By  assuming  a constant  density  for  the 
particles  measured  within  the  analytical  size  range  (4.0  to  10.6  4>),  the 
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volume  percent  data  could  then  be  interpreted  in  terms  of  weight  percent 
per  given  size  range.  Folk's  (1974)  statistical  parameters  for  grain- 
size  distribution  were  used  in  the  final  interpretation  of  these  data 
(see  p.  42  for  specific  equations  used). 

Wash  load  composition 

47.  Suspended  matter  composition  of  the  eroded  bed  load  mater- 
ial was  determined  by  binocular  microscope  (7X)  examination  of  wash  load 
material  retained  on  0.40-pm  Nucleopore  filters.  One  dried  filter  from 
each  experimental  run  was  observed  for  particle  identification  other 
than  silt  and  clay.  Also,  a X-ray  diffraction  analysis  was  performed 
on  two  suspended  sediment  samples  obtained  at  a 30.5-cm  water  depth 
during  experimental  runs  using  Buoy  C type  sediment. 

48.  Sample  preparation  of  the  suspended  matter  for  X-ray 
diffraction  analysis  followed  that  of  Jackson  (1954),  and  X-ray  dif- 
fractogram  interpretation  corresponded  to  that  given  by  Scafe  (1968). 
Removal  of  carbonates  from  the  samples  was  accomplished  by  treatment 
with  5 iN  NaOAc  (sodium  acetate).  Addition  of  30  percent  1^02  (hydrogen 
peroxide)  and  heat  removed  organic  matter  and  Mn02  (manganese  oxide). 

49.  Each  sample  was  then  centrifuged  and  split  into  two 
equivalent  aliquots.  One  was  saturated  with  magnesium  (Mg  ) and  the 
other  with  potassium  (K+)  prior  to  sedimentation  onto  Vycor  (heat- 
resistant)  glass  slides.  After  treatment  of  half  of  the  Mg-saturated 
samples  with  ethylene  glycol.  X-ray  diffraction  analysis  was  performed 
on  glycolated  and  Mq-saturated  samples  at  room  temperature.  The 
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K-saturated  samples  were  analyzed  at  room  temperature,  after  heating  to 


300°C,  and  also  after  heating  to  500°C.  The  X-ray  diffraction  apparatus 
consisted  of  a Phillips  X-ray  diffractometer  set  at  Cu  K - alpha  radia- 
tion 30  kV,  20  mA,  and  a scanning  rate  of  l°/min  spanning  4°  to  30°  20. 


50.  A rectangular  circulating  flume  (Figure  16)  was  designed 
and  constructed  to  model  the  effects  of  unidirectional  turbulent  flow 
upon  sand,  silt,  and  clay  bed  load  material  sampled  from  the  Galveston 
site.  The  flume's  outermost  dimensions  measured  744.2  by  243.8  cm. 

The  actual  water-bearing  channel  measured  61.6-cm  wide  and  71.0-cm 
high.  This  width  minimized  current  boundary  effects  caused  by  the  walls, 
and  this  height  permitted  a water  depth  within  which  parametric  vertical 
profiles  could  be  obtained  during  experimental  runs.  A baffle  system 
was  incorporated  within  two  radii  (inner  radius  30.5  cm,  outer  radius 
91.5  cm)  in  order  to  minimize  helical  streamline  patterns  developed  by 
unrestricted  passage  through  channel  turns.  This  system  consisted  of  the 
placement  of  progressively  larger  radii  sheet  metal  strips  within  the 
channel  corners  (Figure  16). 

51.  The  flume  was  constructed  of  16-gage  sheet  metal 
bent  to  the  channel  dimensions  previously  cited  and  then  welded  to 
rigid  angle  iron  supports.  The  cross-sectional  view  shown  in  Figure  17 
illustrates  the  trough  shape  and  its  containment  by  arc-welded  angle  iron 
supports.  All  metal  was  spray-painted  with  a red  oxide  chlorinated- 
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Schematic  illustration  of  reel 


rubber  primer  and  finally  coated  with  a light  blue  chlorinated-rubber 
base  enamel  to  resist  saltwater  corrosion  and  bedload  material  abra- 
sion. Six  1.3-cm-thick  glass  window  panes,  measuring  58.4  by  106.7  cm 
each,  were  installed  on  the  observation  side  and  sealed  into  place 
using  Silastic  732  adhesive/sealant  and  metal  frame  brackets.  These 
windows  permitted  visual  observation  and  photographic  documentation  of 
dredged  material  discharges,  incipient  bedload  motion,  and  bedform 
characteristics. 

52.  Two  centrifugal  pump  systen s were  installed  to  initiate 
and  control  low-turbulence  flow  currents  wit^‘  the  water  column.  Figure 
17  illustrates  these  systems  in  cross-sectional  view.  The  I ,0-hp  pump 
induced  a maximum  current  speed  of  31  cm/sec  within  the  observation 
channel.  Later  in  the  study,  an  additional  2.0-hp  pump  and  filter 
system  was  installed  to  boost  the  flow  speed  to  a 60-cm/sec  maximum. 

(1 

Gate  valves  installed  on  the  pump  discharge  lines  were  used  to  regulate 
outflow  pressure  and  thus  control  current  speed.  All  piping  was  con- 
structed of  PVC  tubing  to  keep  metal  corrosives  from  entering  the 

system. 

' 

Bed  load  material  introduction 

53.  Approximately  68  kg  of  homogenized  sediment  was  intro- 
duced into  the  flume  using  a 27.9-  by  27.9-  by  45.7-cm  release  box 
suspended  immediately  above  the  air-water  interface  at  the  center  of 
the  observation  channel.  Each  discharge  occurred  in  still  water 
except  one  Block  27  sample  released  during  a 30-cm/sec  current  speed. 
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Photographs  and  stopwatch  timing  of  several  discharges  documented 
turbidity  plume  travel  speed  along  the  bottom  of  the  observation 
cfhannel,  The  plume's  suspended  sediment  concentration  was  monitored 
during  one  Buoy  D sample  release.  Suspended  sediment  dispersed  from 
each  discharge  was  permitted  to  settle  out  for  24  to  36  hr  prior  to  the 
initiation  of  unidirectional  currents  in  the  flume. 

Suspended  sediment  sampling 

54.  Equipment  used  for  collecting  water  samples  during  experi- 
mental runs  consisted  of  0.5-4  sampling  bottles  as  shown  in  figure  18. 
These  bottles  enabled  point  sampling  of  the  ambient  current  with  little 
disturbance  to  the  flow.  Sampling  of  turbidity  plumes  generated  by 
sediment  ejected  into  still  water  was  accomplished  by  installation  of 
plastic  siphon  tubing  on  a depth-ca I i brated  ruler.  Either  method  permit- 
ted withdrawal  of  0.5-4  water  samples  in  60  sec  or  less.  Fach  sample 
collected  was  stored  in  0.5-4  Nalge  plastic  bottles  for  TSM  and  grain- 
size  distribution  analysis. 

55.  Sampling  bottles  were  positioned  at  water  depths  of  15.2, 
30,5  and  45.7  cm  in  the  channel  center  at  the  downstream  sections  of  the 
observation  side.  The  30.5-cm  sampling  depth  was  selected  as  the  primary 
sampling  point  since  mean  flow  conditions  existed  at  this  depth.  Thus, 
plots  of  TSM  concentration  versus  time  were  based  upon  samples  collected 
here. 


56.  A total  of  32  experimental  runs  were  performed  on  the  four 
sediment  mixtures.  Four  of  the  runs  had  a duration  of  longer  than  24  hr. 
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Suspended  sediment  samples  were  taken  intermittently  during  each  run. 

The  time  in  minutes  since  the  previous  current  velocity  increase  or 
decrease  was  recorded  for  every  water  sample  obtained  in  order  to 
accurately  assess  the  changes  in  suspended  sediment  concentration  over 
time.  Appendix  A lists  the  TSM  concentration,  time  after  velocity 
change,  and  location  of  each  water  sample  retrieved  during  experi- 
mental runs. 

Current  speed  measurements 

57.  A Price-Pygmy  water  current  meter  was  used  for  point 
measurements  of  the  unidirectional  current  velocities  for  Buoy  C and 
Block  15  experimental  runs.  Later  in  the  study,  a Model  651  Velmeter 
electromagnetic  current  meter  replaced  the  Price-Pygmy  meter  for  flow 
speed  measurements  during  Block  27  and  Buoy  D experimental  runs.  The 
instrument's  0,1-sec  response  time  and  1.0-mm/sec  sensitivity  permitted 
ample  vertical  measurements  to  be  taken  within  the  water  column.  The 
sensor  consisted  of  a streamlined  probe  measuring  27.94  cm  long  with 
sensing  electrodes  located  7.62  cm  from  the  bottom. 

58.  Vertical  point  measurements  of  the  flow  speed  were 
obtained  during  each  experimental  run  in  order  to  calculate  the  approxi- 
mate bed  shear  (shear  stress  acting  upon  the  sediment-water  interface) 
and  mean  current  speed.  Profiles  taken  with  the  Price-Pygmy  current 
meter  consisted  of  point  measurements  obtained  at  6.1,  15.2,  30.5,  45.7, 
and  57,9  cm  above  the  bed.  Velocity  profiles  obtained  with  the  electro- 
magnetic current  meter  consisted  of  point  readings  measured  at  the 
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following  elevations  above  the  bed:  9.1,  12.2,  15.2,  18.3,  21.3,  24.8 
30.5,  42.7,  and  54.9  cm.  Profile  locations  In  the  observation  channel 
were  recorded  at  positions  charted  in  Figure  19. 
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Bed  shear  stress 

59.  Determination  of  the  shear  stress  acting  upon  the 
sediment-water  interface  was  calculated  using  the  Bedload  Function 
for  sediment  transport  in  open  channel  flows  (after  Einstein  1950): 

i 

U = 5.75  log |o  (30.2  _Y_) 

u*  A (13) 

where  (A)  Is  defined  as  the  apparent  bed  roughness  given  by  (ks/x) 
where  (ks)  is  the  roughness  of  the  bed  and  (x)  is  a corrective  para- 
meter. The  friction  velocity  (u#)  was  determined  by  regression 
analysis  such  that: 

u„  = m/5.75  (14) 

] 

where  (m)  is  defined  as  the  slope  of  plot  of  log  (Y)  versus  (U). 

The  bed  shear  was  computed  from  the  following  eguation: 

tq  = pu*2  (15) 

2 

where  (t  ) is  measured  in  dynes  per  cm  and  (u»)  is  measured  in  cm 
o * 

per  centimeter  per  second.  Computed  values  of  the  bed  shear  stress 
for  each  vertical  profile  are  listed  in  Appendix  B. 

Sal i n I ty  control 

60.  Artificial  seawater  was  made  by  adding  21  g sodium 
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chloride  (NaCL)  and  4 g magnesium  sulphate  (MgSO^)  per  liter  of  tap 
water  to  reproduce  the  near-bottom  salinity  of  28  °/oo  that  exists 
in  the  study  area  (after  Cool  1976),  A salinity  titration  technique 
was  employed  to  monitor  the  salinity  prior  to  experimental  runs. 
Reproduction  of  the  bottom  water  salinity  was  necessary  to  duplicate 
the  degree  of  clay  mineral  flocculation  In  situ. 

i 
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TART  IV:  RESULTS 

Geotechnical  Properties  of  Bed  load 

Grain-size  distribution 

61.  Dispersed  grain-size  analyses  performed  on  sediment 
mixtures  sampled  from  Buoy  C,  Buoy  D,  Block  27,  and  Block  15  stations 
(see  Figure  3)  Indicate  that  all  four  types  had  a median  grain  size 

In  the  silt  range,  were  very  to  extremely  poorly  sorted,  and  were  near 
symmetrical  to  finely  skewed.  The  statistical  measures  used  to  describe 
the  grain-size  distributions  are  given  by  Folk  (1974).  Cumulative 
curves  showing  the  percentage  distribution  of  particle  sizes  plotted  on 
a probability  axis  and  each  phi  size  class  plotted  on  an  arithmetic  axis 
served  as  the  basis  for  these  determinations.  These  curves  are  shown 
in  Figures  20  through  23. 

62.  The  two  control  areas.  Block  15  and  Block  27,  exhibited 
textural  distributions  correlative  with  their  proximity  to  the  Galveston 
barrier  island  system.  The  50.5  percent  sand  content  of  Block  15  mater- 
ial Is  attributable  to  the  proximity  of  sampling  areas  (4.8  km)  to  the 
Galveston  Bay  entrance  channel,  where  Coulthard  (1976)  documented  the 
occurrence  of  a major  sand  source.  This  sand  presumably  is  the  result 
of  high-energy  tidal  currents  within  and  adjacent  to  the  entrance 
channel.  Hall  (1976)  suggested  that  this  control  block  represents  an 
upper  shoreface  sediment  of  a typical  barrier  Island  sequence.  Rlock 

27  material,  obtained  approximately  II  km  southwest  of  the  entrance 
channel,  presumably  represented  a lower  shoreface  zone  sediment  where 
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Figure  20.  Rlock  15  sediment  dispersed  grain-size  distribution 
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Figure  21.  Block  27  sediment  dispersed  qrain-size  distribution 
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Figure  23.  Buoy  D sediment  dispersed  grain-size  distribution 
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silts  and  clays  predominate. 

63.  Figures  22  and  23  illustrate  the  grain-size  distribution 
tor  dredged  material  sampled  at  Buoy  C and  Buoy  D.  The  extremely  poorly 
sorted  nature  (^|  = 4.96)  and  bimodal  distribution  of  Buoy  D material  are 
probably  due  to  the  high  proportions  of  Beaumont  Clay  Formation  sediment 
observed  in  the  samples.  The  high  sand  percentage  (63.7  percent)  ot 
Buoy  C sediment  is  due  largely  to  the  amount  of  shell  debris  (^CO  per- 
cent) observed  in  sand  traction  samples. 

64.  A nondispersed  grain-size  analysis  was  also  run  on  the 
four  sediment  types.  The  following  tabulation  lists  the  sandrsi 1 t:clay 
percentages  and  the  calculated  grain-size  distribution  parameters  pre- 
sented by  Folk  (1974).  In  comparison  to  the  dispersed  analysis,  the 
nondispersing  technique  resulted  in  a decrease  in  both  the  median 
grain  size  (phi  units)  and  inclusive  graphic  standard  deviation  and  an 
increase  in  the  inclusive  graphic  skewness  for  Buoy  C,  Buoy  D,  and 
Block  27  sediment.  Block  15  material  showed  an  increase  in  all  three 
parameter  values.  Substitution  of  distilled  water  for  Calgon  disper- 
sant permitted  clay  flocculation  during  the  pipette  analysis,  thus 
resulting  in  a shift  of  the  grain-size  distribution  to  lower  phi  sizes 
(larger  grain  diameters).  An  opposite  shift  for  Block  15  sediment  is 
probably  the  result  of  incomplete  homogen i zat Ion  of  the  sample  prior 

to  analysis. 
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Buov  C 

Buov  D 

Block  15 

Block  27 

Sand,  % 

64.9 

64.3 

36.7 

13.6 

Silt,  % 

16.5 

11.9 

30.4 

55.1 

Clay,  % 

18.6 

23.7 

32.9 

31  .3 

Mz  (<f>) 

4.80 

3.77 

6.83 

6.53 

*1 

2.61 

4.14 

4.20 

2.50 

SKj 

0.78 

0.21 

0.60 

0.31 

Physical  properties 

65.  The  following  physical  properties  were  measured  for 

sediment  recovered  from  Buoy  C,  Buoy  D,  Block  15,  and  Block  27:  (I) 

void  ratio,  (2)  water  content,  (3)  bulk  density,  (4)  specific  gravity, 

(5)  remolded  shear  strength,  (6)  organic  content,  and  (7)  Atterberg 
limits.  Values  of  these  physical  properties,  along  with  sand : s 1 1 1 :c I ay 
percentages  determined  from  dispersed  grain-size  analysis,  are  presented 
in  Table  2, 

I 

66.  As  anticipaled,  an  increase  in  void  ratio  and  water  content 
corresponded  to  a decrease  in  bulk  density  for  the  four  sediment  mix- 
tures. With  the  exception  of  the  remolded  shear  strength  values.  Buoy 
C and  Buoy  D dredged  material  exhibited  similar  physical  property  values. 
Figure  24  Illustrates  remolded  shear  strength  (Cu>  measurements  via 
plots  of  torque  versus  rotation  angle  for  the  four  mixtures.  The  high 
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shear  strength  of  Buoy  D sediment  is  probably  attributable  to  the 
highly  consolidated  Beaumont  Clay  Formation  present.  The  variations  in 
the  physical  properties  between  Block  15  and  Block  27  material  are  due 
to  the  proportions  of  clay  present  in  the  samples.  The  45.6-percent 
clay  content  in  Block  27  material  accounts  for  the  larger  values  of 
water  content,  void  ratio,  organic  content,  shear  strength,  and  plasti- 
city  index  than  for  Block  15  sediment  containing  only  24.9  percent  clay. 

67.  Values  of  the  plastic  limit  and  plasticity  index  (valid 
for  the  l.25-4>  fraction)  for  Block  15,  Buoy  C,  and  Buoy  0 sediment  are 
defined  by  the  Unified  Soil  Classification  System  (USCS)  (USWES  1955) 
as  "inorganic  fine  sandy  or  silty  soils,  type  MH."  This  is  expected 
since  wet-sieving  the  samples  through  a U.S.  Standard  No.  40  (0.42-mm) 
sieve  removed  particles  larger  than  medium  sand  while  allowing  fine  sand 
and  silt,  if  present,  to  pass  through.  The  values  of  the  plastic  limit 
(63.5)  and  plasticity  index  (38)  for  Block  27  sediment  correspond  to  a 
USCS  classi f ication  as  an  "inorganic  clay  of  high  plasticity,  type  CH." 
This  Is  attributable  to  the  high  clay  content  (45,6  percent)  present. 

68.  Values  of  the  organic  content  and  specific  gravity  ranged 
from  1.84  to  4.14  percent  and  2.67  to  2.71,  respectively.  Compositional 
differences  between  sediments  can  not  be  discerned  on  the  basis  of  these 
low  variations. 

I 

Flume  Experiments 

Bed  I oad  Introduction 

69.  The  discharge  of  each  test  sediment  into  the  water  column 
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resulted  In  simultaneous  turbidity  plume  generation  and  mounding  of 
sediment  on  the  flume  bottom.  The  mounds  exhibited  a hummocky  relief 
up  to  20  cm  high  and  were  spread  about  1.5  meters  upstream  and  down- 
stream from  the  discharge  point.  All  discharges  occurred  in  still 
water  except  for  one  Block  27  material  discharge  released  into  a 
30-cm/sec  flow.  The  resultant  (Block  27)  mound  formed  an  elliptical 
configuration  approximately  1.2  m downstream  of  the  ejection  point. 
Figure  25  illustrates  the  turbidity  plume  generated  by  a 34-kg  Buoy  D 
sediment  discharge.  Stopwatch  timing  of  the  plume  travel  speed 
Immediately  after  ejection  indicated  an  approximate  head  speed  of  7.0 
cm/sec.  Also,  two  water  samples,  taken  during  the  Initial  passage  of 
the  plume  1.2  m downstream  of  the  discharge  site,  were  determined  to 
have  an  average  TSM  concentration  of  50  mg/I  at  a 45-cm  water  depth. 

Erosional  runs 

70.  Four  erosional  runs  were  made  on  both  Buoy  C and  Block 
15  mounds  (and  associated  turbidity  plume  deposits)  using  the  1.0-hp 
centrifugal  pump  system  to  generate  currents  up  to  31  cm/sec.  Using 
an  additional  2.0-hp  pump  system  to  boost  the  flow  speed  to  60  cm/sec, 
n|ne  erosion  runs  were  made  on  both  Block  27  and  Buoy  D sediment.  In 
addition  to  scouring  the  bed  (herein  termed  dense  bed)  that  resulted 
from  the  discharges,  a second  bed  type  was  eroded  during  Block  27  and 
Buoy  D erosional  runs.  This  bed  (herein  called  floe  bed)  consisted 
primarily  of  a flocculated  loose  bed  formed  by  deposition  of  suspended 
matter  scoured  from  the  dense  bed.  A listing  of  all  the  experimental 
runs  In  their  sequence  of  occurrence  Is  contained  In  Appendix  A. 
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Buoy  D Sediment  Discharge 


numbered  in  tenths  of  a foot) 


71.  Dense  bed — sediment  transport.  Determination  of  the  mode 


of  sediment  transport  for  each  of  the  four  sediment  types  was  made  on 
the  basis  of  monitoring  the  TSM  concentration  over  time  per  given  flow 
rate.  Plots  of  the  log  of  the  TSM  concentration  versus  time  per  given 
flow  rate  for  each  experimental  run  are  shown  in  Appendix  C.  For  the 
purpose  of  comparing  the  erodibility  of  the  four  sediment  mixtures, 
selected  log  TSM  versus  time  curves  were  combined.  Figure  26  illus- 
trates such  a graph  for  Buoy  C and  Buoy  D dredged  sediment  erosional 
runs.  As  shown,  the  TSM  concentrations  increased  logarithmically  during 
the  initial  few  hours  following  each  current  speed  increase  and  then 
assumed  a linear  slope  for  the  duration  of  the  experimental  run  until 
the  next  current  speed  change  was  implemented.  The  initial  rapid 
increase  in  the  TSM  concentrations  may  be  attributable  to  resuspension 
of  previously  deposited  turbidity  plume  sediment  in  the  flume  channel, 
in  addition  to  scouring  of  material  from  the  mound.  The  TSM  concen- 
tration decrease  exhibited  by  the  16-cm/sec  Buoy  D run  indicates  that 
some  deposition  of  suspended  matter  occurred.  Comparison  of  Buoy  C 
and  Buoy  D erosional  runs  (Figure  26)  illustrates  the  similarity 
between  the  average  TSM  concentration  held  in  suspension  by  a given 
flow  rate.  This  result  is  not  surprising  since  the  two  dredged 
material  sediments  exhibited  similar  geotechnical  properties  (Table  2). 

72.  Figure  27  illustrates  a comparison  of  Block  15  and  Block 
27  dense  bed  erosional  runs.  The  log  TSM  versus  time  curves  display 
a pattern  similar  to  that  of  the  Buoy  C and  Buoy  D scour  run  curves 
(Figure  26).  Variations  in  the  TSM  concentration  values  for  near- 
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runs 


BLOCK  77 


equivalent  flow  rates  were  caused  by  two  factors.  First,  a different 
pump  outlet  system  was  used  for  Block  15  experimental  runs.  In  lieu 
of  the  open-pipe  exhaust  port  Illustrated  in  Figure  17,  a set  of  12 
0.32-cm-diam  nozzles  was  attached  to  the  1.0-hp  pump  exhaust  pipe 
prior  to  Block  15  experimental  runs.  As  revealed  by  the  velocity 
profiles  and  calculated  bed  shears  listed  in  Appendix  B,  this  system 
apparently  Imparted  a steeper  velocity  gradient  to  the  channel  flow 
that  resulted  in  larger  bed  shears  than  those  produced  by  the  simple 
open-pipe  outlet.  This  explains  the  higher  TSM  concentration  level 
exhibited  by  the  3 1 -cm/sec  Block  15  run  since  the  calculated  bed 
shear  (2.53  dynes/cm.2)  of  this  run  is  greater  by  a factor  of  two  than 
the  bed  shear  (1.12  dynes/cm^)  for  the  30-cm/sec  Block  27  run.  How- 
ever, this  explanation  does  not  account  for  the  observed  TSM  concen- 
tration levels  and  calculated  bed  shears  for  the  23-cm/sec  flow  rate. 
The  second  factor  that  may  have  influenced  the  TSM  concentration 
variations  per  given  flow  speed  is  the  effect  caused  by  the  different 
physical  properties  of  the  two  sediment  types  (Table  2).  Further 
experimentation  using  a singular  pump  system  and  smaller  current 
velocity  increments  between  experimental  runs  is  needed  to  explain 
such  TSM  concentration  variations. 

73.  A distinct  textural  coarsening  of  the  bed  surface  was 
observed  to  form  during  current  flows  exceeding  16  cm/sec.  Apparently, 
winnowing  of  fine  silts  and  clays  from  the  bed  resulted  in  an  armoring 
of  the  bed  surface  by  coarser  silt,  sand,  and  shell  material  that 
remained  Intact.  The  Increased  resistance  to  erosion  caused  by  this 
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grain-size  increase  at  the  sediment-water  interface  probably  reduced 
the  availability  of  the  finer  bedload  fraction  for  resuspension.  This 
may  partially  account  for  the  reduction  in  the  rate  of  TSM  concentration 
increase  as  shown  by  the  linear  segments  of  the  scour  curves.  Quanti- 
tative information  on  physical  changes  of  the  surface  texture  during 
scour  Is  needed  to  verify  this  hypothesis. 


74.  Distinct  textural  interfaces  were  observed  to  form  against 
the  center  glass  pane  in  the  flume.  Figure  28a  illustrates  the  sediment- 
water  interface  16  hr  after  discharge  of  Buoy  D sediment.  The  light-gray 
sediment  on  the  bed  surface  represents  an  oxidized  layer  composed  of 
flocculated  clays  deposited  from  turbidity  plume  suspensions.  Figure 
28b  illustrates  the  same  bed  surface  after  2 weeks  of  experimental  runs. 
The  thin  layer  depicted  as  interface  A formed  as  a result  of  flocculant 
clays  deposited  at  the  conclusion  of  the  last  erosion  run.  Interface  B 
was  the  remnant  oxidized  layer  shown  In  Figure  28a.  A lag  deposit 
comprised  mostly  of  coarse  silt,  sand,  and  shell  hash  was  exhibited 
between  these  two  interfaces.  This  deposit  formed  during  erosional  runs 
exceeding  about  23  cm/ sec.  The  nonuniformity  of  this  deposit  illus- 
trates the  large  textural  variations  that  can  form  within  the  sediment 
column  as  a result  of  scour  over  a hydrodynamlcal ly  rough  surface. 

75.  Dense  bed — critical  shear  stress  and  erosion  velocity. 

A plot  of  the  bed  shear  stress  In  dynes  per  cirr  versus  the  average  TSM 
concentration  in  mg  per  fc  for  dense  bed  erosional  runs  is  illustrated 
in  Figure  29.  The  average  TSM  concentration  per  erosion  run  was  ascer- 
tained by  calculating  the  mean  TSM  concentration  of  the  linear  slope  of 
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Sediment-water  interface  lo  in*  after  discharge 


Interface  A 


Interface  B 


Sediment  interfaces  fol lowing  J weeks  of 
experimental  runs 


Sediment-water  and  sediment  interfaces  following 
discharge  of  and  experiments  with  Ruov  P sediment 
(vertical  scale  is  numbered  in  tenths  of  a foot ) 
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each  log  TSM  versus  time  curve.  Bed  shear  stresses  were  determined  from 
the  velocity  profiles  using  Equations  13,  14,  and  15.  Most  of  these 
profiles  were  measured  at  or  near  the  foreslope  of  each  mound.  Appendix 
B lists  the  velocity  profiles,  resultant  bed  shears,  and  profile  loca- 
tions obtained  during  experimental  runs. 

76.  Figure  29  illustrates  that  for  the  sediment  types  tested 
the  suspended  sediment  concentration  increased  rapidly  for  bed  shears 
exceeding  approximately  1.0  dyne  /cm  . This  bed  shear  value  then 
represents  the  critical  shear  stress  (Tocr)  necessary  to  initiate  rapid 
erosion  of  the  dense  bed.  Using  Einstein's  Bedload  Function  for  open 
channel  flow  conditions  (see  paragraph  59),  a 1.0-dyne  /cm^  bed  shear 
corresponds  to  a critical  erosion  velocity  of  24  cm/sec  at  30  cm  above 
the  bed.  This  velocity  agrees  with  visual  observation  of  mass  bed  load 
transport  for  current  speeds  exceeding  23  cm/sec  measured  at  30.5  cm 
above  the  bed. 

77.  The  similarity  of  the  critical  shear  stress  value  for 
the  four  sediment  types  suggests  that  the  degree  of  cohesiveness 
between  sedimentary  particles  Is  nearly  equivalent  for  these  four 
sand,  silt,  and  clay  mixtures.  This  is  true  assuming  that  the  resis- 
tance to  erosion  is  primarily  dependent  upon  the  degree  of  cohesive 
force  acting  between  sedimentary  particles. 

78.  Floe  bed.  In  order  to  examine  the  erosional  mode  of 
sediment  transport  on  a loosely  deposited  flocculated  clay  bed,  three 
scour  runs  were  performed  on  both  Block  27  and  Buoy  D sediment  deposited 
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from  suspension.  The  first  of  these  runs  was  initiated  after  a settling 
period  that  followed  the  last  deposition  run.  The  14-cm/sec  Block  27 
floe  bed  erosion  run  was  initiated  after  a 36-hr  settling  period,  whereas 
the  16-cm/sec  Buoy  D floe  bed  erosion  run  began  following  a 48-hr  set- 
tling period.  In  both  cases,  the  channel  bottom  and  sediment  mound  were 
uniformly  blanketed  with  a thin  ( 0.32-cm-thick)  layer  composed  of  floc- 
culant  clays  and  very  fine  silts. 

79.  The  scour  curves  for  these  erosional  runs  are  exhibited 
in  Figures  30  and  31.  As  shown,  the  TSM  concentration  tended  to  in- 
crease rapidly  to  some  asymptotic  level.  A slight  decrease  in  the  TSM 
concentration  with  time  was  evident  for  Block  27  Run  No.  9 and  Buoy  0 
Run  Nos.  9 and  10,  which  Indicated  that  deposition  occurred.  Much  of 
the  material  that  was  resuspended  probably  originated  as  a floe  bed 
in  the  return  channel  of  the  flume.  Since  the  current  propulsion  system 
was  located  here,  much  higher  bed  shears  were  produced  than  in  the 
observation  channel.  This  was  substantiated  by  the  higher  TSM  concen- 
tration levels  attained  by  the  floe  bed  runs  than  for  the  dense  bed 
runs  made  at  similar  flow  rates.  Unfortunately,  the  lack  of  point 
velocity  measurements  within  the  return  channel  prevented  an  accurate 
estimation  of  the  critical  bed  shear  required  to  induce  floe  bed  load 
erosion.  However,  the  data  collected  suggest  that  the  pattern  of  floe 
bed  erosion  was  similar  to  that  exhibited  by  the  dense  bed  erosion 
runs;  that  is,  the  TSM  concentration  of  the  fluid  appears  to  be  directly 
dependent  upon  the  magnitude  of  the  bed  shear  stress  acting  on  the 
sediment-water  Interface. 
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80.  Erosion  rates.  Due  to  the  hydrodynamica I I y rough  surface 


produced  by  mounding  of  sediment  on  the  flume  bottom,  an  accurate  esti- 
mate of  the  erosion  rate  could  not  be  ascertained;  that  is,  the  nonuni- 
formity of  the  bed  surface  (Figure  28b)  prevented  accurate  determination 
of  the  actual  surface  area  subjected  to  scour.  In  spite  of  this,  an 
approximation  of  the  relative  scour  rate  was  gained  from  an  evaluation 
of  the  rate  of  TSM  concentration  increase  as  determined  from  the  linear 
slope  of  the  erosion  curves.  This  slope  was  given  by  a visual  regres- 
sion line  drawn  through  the  TSM  concentration  values  measured  between 
4 to  24  hours  after  each  velocity  change. 

81.  Table  3 lists  the  rate  of  change  of  the  TSM  concentration 
per  average  flow  speed  for  both  dense  and  floe  bed  eroslonal  runs.  The 
maximum  observed  rate  increase  given  by  the  linear  slope  was  16.47 
mg/l/hr  for  an  average  sustained  flow  speed  of  57  cm/sec  on  Block.  27 
sediment.  In  comparison  to  the  values  listed  in  Table  3,  this  TSM 
concentration  rate  increase  seems  high.  However,  the  naturally  high 
clay  content  (45.6  percent)  of  this  lower  shoreface  sediment  would  be 
expected  to  result  in  higher  suspended  clay  concentrations  during  ero- 
sion than  for  a bedload  sediment  containing  less  clay  material. 

82.  For  flow  speeds  ranging  from  10  to  44  cm/sec,  observed 
TSM  concentration  rate  changes  ranged  from  -0.36  to  2.17  mg/l/hr. 

Such  low  rates  indicated  that  the  relative  scour  rate  under  low  to 
medium  flow  regime  conditions  was  small.  Observations  of  turbidity 
plume  formation,  coupled  with  high  rates  of  TSM  concentration  increase 
during  the  first  few  hours  of  scour,  indicated  that  the  bulk  of  clay 
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Table  3 
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Rate  of  Change  of  TSM  Concentration  For  Erosional  Runs 


Rate  of  Change  of 

Type 

Average  Current 

TSM  Concentration 

Buoy  Of  Bed 

Run  No. 

Velocity  (cm/sec) 

mg/i/hr 

Buoy  C Dense 

1 

10 

0.09 

2 

15 

0.25 

3 

23 

0.51 

4 

32 

1 .02 

Dense 

1 

16 

-0.36 

2 

23 

0.01 

3 

32 

0.24 

4 

39 

0.40 

Buoy  D 

5 

44 

0.50 

6 

60 

1 .42 

Floe 

9 

16 

-0.23 

10 

24 

0.23 

1 1 

31 

1 .80 

Block  15  Dense 

2 

16 

0.02 

3 

23 

-0.04 

4 

31 

0.74 

Dense 

1 

12 

-0.10 

2 

23 

0. 1 1 

3 

30 

0.14 

4 

38 

1 .06 

Block  27 

5 

44 

2.17 

6 

57 

16.47 

Floe 

9 

14 

-0.06 

10 

22 

1 .10 

1 1 

30 

1 .73 
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and  silt  resuspension  occurred  during  the  first  few  hours  immediately 
following  a dredged  material  discharge.  Then  scour  rates  were  reduced 
to  small  values  as  textural  coarsening  of  the  bed  surface  occurred  and 
perhaps  also  as  the  sediment  consolidated  under  the  effect  of  gravity. 

Deposit  Iona  I runs 


83.  Experiments  were  performed  to  determine  the  approximate 
current  velocity  and  associated  bed  shear  at  which  significant  deposi- 
tion of  the  suspended  sediment  occurred.  Figure  32  illustrates  the 
depositional  runs  made  using  each  of  the  four  sediment  mixtures.  The 
deposition  runs  for  Block  15  and  Buoy  C sediments  were  initiated  at 
the  conclusion  of  the  32-cm/sec  erosional  runs,  whereas  Block  27  and 
Buoy  D deposition  runs  began  immediately  following  the  57-  and  60-cm/sec 
erosion  runs,  respectively.  In  this  manner,  the  suspended  sediment 
concentration  that  resulted  from  scour  of  the  dense  bed  was  monitored 
over  time  at  a 30.5-cm  depth. 

84.  As  shown  in  Figure  32,  significant  deposition  of  the 
suspended  sediment  occurred  during  sustained  flow  speeds  less  than 
about  10  cm/sec  and  initial  flow  concentrations  exceeding  approximately 
100  mg/4.  Deposition  curves  for  the  6-cm/sec  Block  27  run  and  4-cm/sec 
Buoy  D run  Indicate  that  the  suspended  sediment  concentration  decreased 
logarithmically  for  the  24-hr  period  of  observation . For  current 
velocities  exceeding  10  cm/sec,  the  TSM  concentration  remained  stable 
indicating  that  suspended  silts  and  clays  are  apt  to  remain  in  suspen- 
sion indefinitely.  An  average  flow  speed  of  10  cm/sec,  30  cm  above  the 
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bed,  corresponds  to  a calculated  bed  shear  of  approximately  0.20 
2 

dyne  /cm  as  determined  from  Einstein's  Bed  load  function  for  open 
channel  flow  conditions  (paragraph  59).  This  value  should  then  repre- 
sent the  lowest  bed  shear  which  inhibits  significant  deposition  of 
suspended  sediment.  This  value  agrees  well  with  the  lowest  bed  shear 
(0.23  dyne  /cm  ) above  which  rapid  suspended  sediment  deposition  did 
not  occur  in  the  flume  (see  experimental  run  No.  6,  Buoy  C,  in  Appendix 
B). 

Geotechnical  properties  of  wash  load 

85.  Grain-size  distribution.  The  grain-size  distributions  of 
the  sediment  suspended  in  water  samples  obtained  during  experimental 
runs  were  determined  over  a total  size  range  of  4 . 0<}>  to  1 0 . 6<t>  (coarse 
silt  to  clay).  Cumulative  curves  plotted  from  the  measured  cumulative 
weight  percent  per  size  increment  served  as  the  basis  for  determination 
of  the  graphic  mean.  The  resultant  values  for  each  analysis  are  tabu- 
lated in  Appendix  D,  in  addition  to  sample  location,  mean  flow  speed, 
time  of  sampling,  and  experimental  run  designation.  The  grain-size 
distribution  results  indicate  a range  of  the  graphic  means  from  6.5  to 
9 . 7<J>  (fine  silt  to  clay)  and  an  average  clay:silt  ratio  of  3.9:1. 

86.  Cumulative  curves  exhibiting  the  grain-size  distribution 
of  three  water  samples  obtained  during  Block  15,  dense  bed  31-cm/sec 
erosional  run  are  illustrated  in  Figure  33.  A slight  decrease  in  the 
silt  fraction  with  a simultaneous  increase  in  the  clay  fraction  is 
indicated  for  samples  obtained  during  increasing  time  after  velocity 

change.  A similar  pattern  Is  depicted  in  Figure  34  for  three  samples 
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Figure  33.  Grain-size  distributions  of  three  washload  samples  collected 
during  Block  15,  dense  bed  31-cm/sec  erosion  run 
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CUMULATIVE  WEIGHT  PERCENT 
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Figure  34.  Grain-size  distributions  of  three  washload  samples  collected 
during  Buoy  D,  dense  bed  23-cm/sec  erosion  run 
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collected  during  Buoy  D,  dense  bed  23-cm/sec  eroslonal  run.  Both  sets 
of  samples  were  recovered  30  cm  above  the  bed.  The  decrease  In  the  silt 
content  for  these  sustained  flow  rates  may  have  resulted  from  the  deposi- 
tion of  silt  that  originated  as  resuspended  turbidity  plume  deposits  In 
the  return  channel.  The  higher  shearing  rates  produced  within  the 
return  channel  by  the  current  propulsion  system  probably  caused  rapid 
resuspension  of  the  silt  and  clay  fraction  plume  material  during  each 
Initial  current  speed  increase.  Then,  lower  bed  shears  In  the  observa- 
tion channel,  which  resulted  from  frictional  drag  along  the  channel 
surfaces,  allowed  deposition  of  those  silt-size  fractions  that  exceeded 
the  stream  competency.  However,  further  experimentation  is  desired  to 
determine  the  magnitude  of  the  shearing  rate  In  the  return  channel  per 
given  flow  speed,  and  the  correspondent  grain  size  of  the  resuspended 
mater i a I . 

87.  Large  variations  in  the  washload  grain-size  distribution 
were  noted  for  water  samples  collected  at  different  flow  rates.  Figure 
35  Illustrates  the  textural  and  compositional  differences  between  samples 
obtained  during  Block  27,  dense  bed  44-,  57-,  and  30-cm/sec  experimental 
runs.  These  samples  were  also  collected  30  cm  above  the  bottom.  As 
revealed,  a 30  to  77  percent  increase  In  the  silt  content  occurred  after 
the  mean  flow  speed  was  increased  from  44  to  57  cm/s ec.  Next,  a 77  to 
21  percent  silt  content  decrease  took  place  as  the  flow  rate  was  reduced 
to  30  cm/ sec.  Another  example  displaying  an  increase  in  the  silt  content 
with  increasing  current  speed  Is  shown  In  Figure  36  for  Block  27,  floe 
bed  erosion  runs.  These  silt-content  variations  illustrate  the  depen- 
dence of  the  washload  grain-size  distribution  upon  the  competency  of  a 
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Figure  35.  Grain-size  distributions  of  three  washload  samples  collected 
during  Block  27,  dense  bed  44-,  57-,  and  30-cm/sec  erosional 
runs 
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Figure  36.  Grain-size  distributions  of  three  wash  load  samples  collected 
during  Block  27,  floe  bed  I4-,  22-,  and  30-cm/sec  eroslonal 
runs 


97 


given  flow.  As  revealed  by  these  examples,  an  Increase  or  decrease  In 
the  blmodallty  of  the  groin-size  distribution  occurs  with  an  Increase 
or  decrease  In  the  flow  rate. 


88.  Composition,  Binocular  microscope  examination  of  sus- 
pended sediment  retained  upon  0.40-pm  Nucleopore  filters  Indicated  that 
fine  silt  and  clay  comprised  the  bulk  of  the  suspended  matter.  Occa- 
sional organic  fragments  ( foraml n I fora  I tests)  and  siliceous  grains 
were  noticed  on  a few  of  the  filters,  but  their  abundance  did  not  exceed 
I or  2 percent  of  the  total  composition.  The  most  prevalent  component 
present  on  the  filters  was  clay.  As  stated  in  the  previous  section, 
washload  grain-size  analyses  Indicated  a net  c I ay : s 1 1 1 ratio  of  3.9:1, 
which  confirms  this  observation. 

89,  An  X-ray  diffraction  analysis  for  clay  mineral  deter- 
mination was  performed  on  two  suspended  sediment  samples  recovered 
during  two  Buoy  C experimental  runs.  Figures  37  and  38  Illustrate 
X-ray  dl f f ractograms  of  suspended  matter  sampled  from  a 3 1 -cm/sec  ero- 
sion and  23-cm/sec  deposition  run,  respectively.  Clay  minerals  were 

Identified  by  their  first-order  angstrom  spaclngs.  Montmor I I Ion  I te 

o 

Is  characterized  by  a 17.6-A  peak  on  an  ethy lene-glycolated  and 

Mg-saturatod  sample.  Also,  It  is  Identified  by  collapse  to  10-^  with 

heating  of  a K-saturated  sample.  Both  chlorite  and  vermlcullte  display 

a 14.2-X  peak.  Illlte  exhibits  a 1 0— A Mg-glycolated  peak  that  increases 

In  height  with  K-saturatlon  plus  heat,  A chlorl te-kaol Inlte  peak 
o 

occurs  at  7.2  A on  the  Mg-saturated  and  glycolated  dl tf ractogram. 
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90.  Estimates  of  the  relative  percentages  of  the  clay  minerals 
were  calculated  using  the  peak  height  method  described  by  Scale  (I9t>8). 
The  I I I i te  peak  (10  X)  was  utilized  as  the  internal  standard  In  the  com- 
putations. The  relative  percentage  of  the  clay  minerals  tor  the  less 
than  J-um  traction  of  the  two  samples  are  presented  In  the  following 
tabulation: 


Compos  1 1 1 on  t 


Run 

No. 

lllite 

Kaol Inite- 
Chlor i te 

Montmor i 1 Ion i te 

Ch lor i te- 
Vermicu 1 i te 

4 

59 

25 

15 

5 

5 

58 

21 

17 

4 

As  shown,  lllite  exhibited  the  highest  percent  composition  in  the  samples 
followed  by  kaol i n I te-ch lor i te,  montmor 1 1 Ion i te,  and  ch lor i te-verml cu I i te 
in  decreasing  order  of  abundance.  Cool  (1976)  reported  high  i I I i te  per- 
centages in  both  suspended  matter  and  bottom  sediment  samples  collected 
from  the  offshore  disposal  site.  He  suggested  that  the  source  of  this 
i 1 1 i te  was  the  Pleistocene  Beaumont  Clay  Formation  that  is  subaeriallv 
exposed  along  portions  of  the  Calveston  Ship  Channel.  Thus,  scour  of 
Buoy  C sediment,  which  presumably  was  dredged  from  the  entrance  channel, 
resulted  In  an  expected  high  lllite  content  In  the  two  wash  load  samples. 

Summary 

91.  Laboratory  flume  experiments  performed  on  the  tour  sedi- 
ment mixtures  Indicated: 


8 

I 

1 

ii 
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Each  of  the  four  sediment  mixtures  eroded  similarly  as 

evidenced  by  the  s i • lar  shape  of  the  scour  curves  and 

2 

a similar  critical  shear  stress  value  (^=1.0  dyne  /cm  ) 
necessary  to  cause  both  massive  bed  load  transport  and 
rapid  TSM  concentration  increases. 

Primary  silt  and  clay  resuspension  occurred  during  the 
discharge  process  and  the  initial  2 to  3 hr  immediately 
following  each  current  speed  increase. 

The  winnowing  of  silts  and  clays  from  the  bed  by  erosion 
resulted  in  a textural  coarsening  of  the  bed  surface. 
Continuous  scour  over  a hydrodynamical ly  rough  bed  pro- 
duced textural ly  diverse  sediment  interfaces  within  the 
sediment  column. 

Rapid  deposition  of  suspended  fine  silts  and  clays 
occurred  at  bed  shears  less  than  approximately  0.2 
dyne  /car  for  initial  TSM  concentrations  exceeding 
about  100  mg/5.. 

Fine  silt  and  clay  comprised  the  bulk  of  suspension 
produced  by  scour.  The  wash  load  grain-size  distribution 
depended  primarily  upon  the  competency  of  a given  flow. 
An  increase  or  decrease  in  the  silt  mode  of  the  distri- 
bution occurred  with  an  Increase  or  decrease  in  the  flow 
rate. 

Illite  was  the  primary  clay  component  produced  by  scour 
of  Buoy  C dredged  material. 


Hydrography  of  Offshore  Galveston 


Tide  regime 

92.  Tide  readinas  obtained  from  the  U.S,  Army  Engineer 
District,  Galveston,  were  recorded  at  the  south  jetty  lighthouse  ot 

the  Galveston  Bay  entrance  channel.  These  readings  consisted  ot 
continuous  strip  chart  records  taken  from  February  1975  through 
I May  1976.  Normal  tides  ranged  from  0.3  to  0.6  m In  height,  and 
were  mixed  (l.e.,  diurnal  and  semidiurnal  in  frequency).  Figure 
39  illustrates  a typical  mixed  tide  cycle  for  the  period  7 through  21 
July  1975. 

93.  Adverse  meteorologica I conditions  were  found  to  alter 

tide  heights.  Hourly  wind  speed  and  direction  data  for  the  time 
period  February  1975  through  May  1976  were  obtained  from  the  Galves- 
ton Weather  Bureau.  Examination  of  both  tide  and  meteorologic  data 
revealed  that  primary  storm  wind  azimuths  affecting  tide  heights 
include  (I)  strong  persistent  east,  southeast,  or  south  winds  blowing 
across  the  Gulf  of  Mexico  for  several  days  and  (2)  strong  persistent 
north  to  northwest  winds  following  the  passage  of  a cold  front.  In 
the  first  instance,  tides  up  to  0.9  m above  normal  were  observed, 
while  in  the  second  case,  tides  up  to  0.6  m nelow  normal  were  recorded. 
Figures  40  and  4|  illustrate  occurrences  of  these  two  phenomena.  It 
Is  presumed  that  strong,  long-duration  southerly  winds  pile  surface 
water  onto  the  coast  resulting  in  a measureable  sea  level  increase, 
whereas  strong  northerly  winds  push  water  seaward  causing  a sea  level 
decrease. 
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PREDICTED  12.5-1 


22  through  25  February  1975 


Current  data 


94.  Numerous  velocity  profiles  and  bottom  current  measure- 
ments were  recorded  to  ascertain  the  circulation  pattern  within  the  study 
area.  Table  I lists  the  type,  location,  date,  and  time  period  of  the 
current  meter  data  recovered.  Data  on  the  tidal  currents  passing  through 
the  entrance  channel  were  discussed  by  Hall  (1976). 

Velocity  prof i les 

95.  Velocity  profiles  obtained  by  Hall  (1976)  at  Buoy  B, 
located  within  the  entrance  channel,  indicated  that  current  speed  and 
direction  are  controlled  by  the  cyclic  passage  of  tidal  currents  through 
the  inlet.  Measured  flood  and  ebb  tide  current  velocities  in  excess  of 
100  cm/sec  were  recorded  during  normal  tide  cycles  within  the  entrance 
channel.  In  contrast,  velocity  profiles  taken  near  Buoys  B,  C,  and  D 
indicate  maximum  current  velocities  of  only  75  cm/sec. 

96.  Vertical  profiles  were  made  at  2-hr  intervals  during  a 
24-hr  anchor  station  at  Buoy  C on  I - 2 May  1976  (during  a diurnal  tide). 
Figures  42  through  44  illustrate  a succession  of  six  profiles  recorded  at 
2000,  2200,  and  2400  hr  I May  and  0200,  0400,  and  0600  hr  2 May.  Three 
important  results  were  gained  from  an  evaluation  of  these  profiles.  First, 
the  speed  and  direction  of  the  currents  in  the  upper  half  of  the  water 
column  are  coincident  with  entrance  channel  tidal  currents.  Second,  a 
northeasterly-directed  bottom  flow  persisted  throughout  the  observation 
period.  Third,  the  near-bottom  sampling  points  In  each  profile  suggest 
either  a shear  zone  between  the  near-bottom  and  upper  water  column  cur- 
rents or  a flow  reversal  near  the  bed. 
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97.  The  1 1 me  of  slack  water  at  the  south  jetty  lighthouse 


(on  the  entrance  channel)  occurred  at  0200  hr,  2 May  I97o,  with  ebb  tide 
preceding  and  flood  tide  following.  Figures  42  and  4'  show  the  southerly 
directed,  near-surface  ebb  flows  prior  to  slack  water  development.  Between 
2400  hr  and  0200  hr,  a sharp  decrease  (from  SO  to  7 cm/sec)  In  near-surface 
flow  speed  occurred  which  corresponded  to  slack  water  during  mean  low  tide 
(Figure  43).  Later,  a north-northeaster  I y t lood-d I reefed  surface  current 
was  initiated  as  Illustrated  in  Figure  44.  Maximum  and  minimum  near-surface 
flow  speeds  varied  with  the  entrance  channel  tidal  range,  as  did  the  sur face 
current  direction.  Also,  during  this  24-hr  observation  period,  sustained 
h-m/sec  southeasterly  winds  occurred.  Such  onshore-d i rected  winds  account 
for  the  clockwise  deflection  of  the  near-surtace  current  direction  with 
depth  due  to  Ekman  drift  (from  Ekman  1 90S). 

90.  In  order  to  compare  current  speed  fluctuations  during 
a tidal  cycle,  a plot  of  flow  velocity  at  various  depths  versus  time  was 
made  using  the  profiles  collected  from  1400  hr  I May  to  1400  hr  2 May  1 97c 
at  Buoy  C.  This  graph  Is  depicted  In  Figure  4S.  Vertical  lines  of  demar- 
cation Indicate  predicted  times  of  slack  water  tor  the  entrance  channel 
as  listed  In  U.S.  Dept,  of  Commerce,  "Tidal  Current  Tables  |9'c>."  As  shown, 
high-speed  surface  ebb  currents  are  coupled  with  low-speed  ebb  flows  lower 
In  the  water  column.  Conversely,  flood-directed  currents  exhibit  greater 
velocities  near  the  bottom  than  at  the  surface.  Also,  the  time  of  flood 
current  commencement  near  the  seabed  precedes  that  In  the  upper  water  col- 
umn. This  suggests  that  a flood-tide  underflow  occurred  near  the  bottom. 
Ludwick  (1974)  also  observed  such  underflows  within  the  tidal  entrance  to 
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00,  lyplcfll  velocity  profiles  obtained  near  Buoy  B and  Buoy 
P are  shown  in  Figures  4o  and  47.  The  profile  recorded  at  Buoy  B at  lOkh 
hr  o April  |07b  exhibits  current  velocities  and  directions  that  are  concur- 
rent with  an  ebb  flow  trow  Calves ton  Bay.  The  Buoy  0 profile  illustrated 
in  Figure  47  Indicates  an  ebb-directed  surface  flow  and  a f I ood-di reefed 
bottom  flow.  This  apparent  flow  reversal  demonstrates  the  spatial  and 
temporal  dissimilarity  between  upper  and  lower  water  column  currents  near 
the  southern  margin  ot  the  study  area.  In  contrast,  profiles  taken  near 
Buoy  B (northern  margin)  suggest  uniformity  between  upper  and  lower  flows. 
This  may  be  explained  by  the  close  proximity  (A.,’  km)  of  the  Buoy  B sta- 
tion to  the  entrance  channel,  where  tide  currents  may  regulate  the  entire 
water  column. 

i 

Offshore  tidal  current  circulation  pattern 

100.  The  tidal  current  regime  exhibited  by  the  Buoy  B,  C,  and 
0 velocity  profiles  Is  consistent  with  the  general  circulation  pattern 
described  by  Hall  (I07p).  Figure  40  Illustrates  the  flow  pattern  developed 
during  flood  tide  at  the  entrance  channel.  The  greatest  current  speeds  and 
subseguent  water  transport  occur  within  the  entrance  channel . Flood  tide 
currents  In  the  study  area  are  typically  oriented  east-northeast.  Figure 
40  Illustrates  the  ebb-flow  circulation  pattern.  In  this  case,  seaward- 
moving  waters  are  deflected  toward  the  southwest  by  prevailing  southeast- 
erly winds  and  the  Coriolis  force  (Hall  |07cO,  Thus,  near-surface  ebb 
tidal  flows  travel  southwest  or  down  coast  parallel  to  Oalveston  Island. 
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Figure  49.  Potary  offshore  flow  pattern  developed  durino  ebb  tide  at  the  Galves 
ton  Bay  entrance  channel  (from  Hall  1976) 


101.  To  evaluate  the  bottom  current  circulation  pattern 


operative  In  the  field  area,  progressive  vector  diagrams  were  prepared 
from  the  continuous  bottom  current  measurements  recorded  at  Ruoy  B and 
Buoy  D.  These  diagrams  were  made  by  computing  a succession  of  time 
integrals  of  horizontal  measurements,  each  starting  from  the  same  Initial 
time.  Hence,  a plot  of  the  resultant  vectors  portrays  a trajectory  with 
spatial  dimensions.  True  water  mass  movement  is  not  actually  illustrated 
by  the  vector  path,  unless  the  currents  lack  spatial  change  over  the  scale 
of  the  path.  Such  spatial  homogeneity  has  not  been  observed  in  the  field 

'1 

area;  thus,  these  diagrams  are  primarily  of  qualitative  value  in  assessing 
net  bottom  water  movement. 

102.  Figure  50  (after  Hall  1976)  illustrates  a progressive 
vector  diagram  plotted  from  15-mln  interval  current  velocity  measurements 
recorded  1.0  m off  the  bottom  at  Buoy  D.  The  period  of  observation  extends 
from  9 October  through  14  November  1975.  Each  plus  sign  (+)  marks  the 
start  of  a 6-hr  period.  Assuming  spatial  homogonolty  of  the  bottom  currents, 

Hall  (1976)  computed  a net  drift  of  260  km  toward  the  southwest. 

103.  Figures  51  through  54  illustrate  progressive  vector 
diagrams  constructed  from  15-mln  interval  bottom  current  readings  recorded 
at  Buoy  B and  Buoy  D during  March  through  April  1976.  The  times  of  obser- 
vation for  these  plots  are  tabulated  In  Table  I.  Also,  calendar  dates  are 
marked  along  the  vector  trajectories  in  each  diagram. 

104.  Two  significant  features  are  illustrated  in  these  diaqrams: 

I IB 


Figure  51.  Progressive  vector  diagram  for  17  March  through 
23  March  1976  at  Buoy  B 
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1200  hr  April  «,l»76 


Figure  53.  Progressive  vector  diagram  for  8 April  through 
30  Apr i I 1976  at  Buoy  B 
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Figure  54.  Progressive  vector  diagram  tor  8 April  through  3 May 
1976  at  Buoy  D 
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(I)  persistent  down-coast  or  offshore-directed  (SW  or  SE)  vector  trajec- 
tories and  (2)  Intermittent  loops.  The  latter  may  represent  the  flpw  rever- 
sals influenced  by  flood  and/or  ebb  tidal  currents,  while  the  former 
signify  the  occurrence  of  onshore-directed  storm  winds.  Examples  of 
these  two  phenomena  fol low. 

105.  The  effect  of  storm-generated  winds  upon  sea  level 
height  is  illustrated  in  Figure  40.  Apparently,  sustained  6.0-  to  10.0- 
m/sec  southeasterly  winds  during  the  5-day  period  14-19  April  1976 
resulted  in  a tide  range  0.6  m above  normal.  During  this  period  a down- 
coast  and  offshore-directed  bottom  flow  persisted  as  revealed  by  the 
time  synchronous  SE-SW  progressive  vector  trajectories  in  Figures  53  and 
54  for  Buoys  B and  D.  These  onshore-directed  storm  winds,  blowing  across 
the  Gulf  of  Mexico  caused  a buildup  of  surface  water  along  the  coast, 
which,  in  turn,  produced  a nearshore  pressure  anomaly  within  the  water 
column.  To  compensate  for  this  positive  anomaly,  mass  continuity  requires 
that  an  opposing  flow  occur  lower  in  the  water  column.  Internal  stratifi- 
cation has  been  demonstrated  by  Murray  (1975)  to  govern  the  location  in  the 
water  column  where  this  opposite  flow  component  exists.  That  a shear  zone 
periodically  exists  near  the  bottom  In  the  disposal  site  Is  evidenced  by 
(I)  observed  flow  reversals  near  the  bed  (Figures  42  through  44);  (2)  con- 
ductivity, salinity,  and  temperature  profiles  recorded  during  May  and  July 
1975  (Cool  1976);  and  (3)  sustained  offshore-directed  bottom  flows  shown 
by  the  SE-SW  vector  trajectories  in  Figures  50  through  54.  Thus,  it  is 
hypothesized  on  the  basis  of  these  data  that  a two- layer  nearshore  circu- 
lation pattern  forms  during  the  passage  of  onshore-directed  storm  winds; 
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that  is,  persistent  southeasterly  winds  exceeding  about  5 m/sec  cause 
onshore  flow  at  the  surface  and  stable  down-coast  to  offshore  currents 
near  the  bottom, 

106.  Another  example  of  this  phenomena  is  Illustrated  by 
the  time-synchronous  wind  and  bottom  current  vectors  plotted  as  stick 
diagrams  In  Figure  55.  The  period  of  observation  extended  from  4 May 
through  8 May  197b.  Doth  wind  and  current  azimuths  are  oriented  toward 
the  direction  of  traverse.  The  bottom  current  vectors  represent  discrete 

* 

2-hr  interval  current  speed  and  direction  (true  north)  measurements 
recorded  at  Ruoy  P and  Duoy  D I m above  the  bottom.  Note  that  a N-NF  to 
S-SW  bottom  current  reversal  occurred  at  both  stations  during  5 May.  On 
this  same  day,  a 0.2-m  positive  tide  anomaly  was  Initiated  as  shown  by 
the  observed  versus  predicted  tide  cycle  diagram  in  Figure  56.  Sustained 
b-  to  iO-m/sec  southeasterly  winds  that  occurred  during  b May  (see  wind 
vector  plot,  Fig.  55)  probably  induced  this  anomaly  via  the  mechanism 
previously  cited.  In  addition,  note  that  the  tide  anomaly  continued  tor 
the  duration  of  the  observation  period  as  did  the  offshore-directed 
bottom  currents.  The  50°  to  90°  offset  between  the  simultaneous  bottom 
current  directions  displayed  at  Buoys  R and  0 may  be  attributable  to 
topographic  effects  created  by  dredged  material  mounding  adjacent  to  the 
buoys  or  perhaps  to  a rotary  bottom  flow  pattern  Indigenous  to  the  area. 
Further  hydrographic  reconnaissance  Is  required  to  explain  such  variations. 
The  Important  point  Is  that  positive  tide  anomalies  produced  by  onshore- 
directed  storm  winds  In  excess  of  about  5 m/sec  generate  relatively 
steady  down-coast  to  offshore-directed  bottom  flows. 
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107.  The  loops  displayed  In  The  vector  trajectories  are 
presumed  Indicative  of  the  influence  of  reversing  tidal  currents  of  the 
entrance  channel.  For  example,  consider  the  bottom  current  speed  and 
direction  diagram  exhibited  in  Figure  bl  for  the  Buoy  0 station.  The 
observation  period  extended  from  ‘>-?o  April  l07o.  Piscrete  '-hr  Interval 
current  speed  and  direction  measurements  were  plotted  from  the  l*i-min 
Interval  data  recorded  by  the  Savonius  current  meter  positioned  I m 
above  the  seabed.  The  time-synchronous  progressive  vector  trajectory  Is 
Illustrated  In  Figure  *’4.  it  is  apparent  that  the  tidal  I y Induced  current 
reversals  depicted  In  Figure  r'7  for  the  per  iods  ll-l\  21-??,  and  7s 
April  I07o  correspond  to  the  SL-SW  trending  loops  exhibited  In  Figure  *'4. 
An  additional  striking  feature  revealed  by  a comparison  of  these  figures 
Is  the  higher  bottom  current  velocities  displayed  during  the  period  of 
southeaster  I y storm  winds  (i.e.,  14-10  April)  than  tor  the  period  where 
reversing  tidal  currents  persisted. 

I OP.  Strong  northerly  winds  also  affect  the  offshore  bottom 
current  regime.  Note  the  large  SI -NN  trending  loop  exhibited  In  I igure 
*>?  tor  the  S-day  period  from  '0  March  to  ' April  ll'7b.  The  observed 
tidal  cycle  is  Illustrated  in  Figure  bp,  and  the  surface  wind  speed  and 
a/imuth  tor  this  period,  i lotted  as  a stick  diagram,  are  shown  in  I Igure 
SO,  wind  azimuths  are  oriented  toward  the  direction  of  travel.  Recon- 
struct ion  of  events  Indicates  that  a strong  northerly  wind  that  occurred 
50-'l  March  effectively  lowered  sea  level  (and  therefore  tidal  range) 
to  a normal  level  by  drawing  surface  water  seaward.  Purl ng  this  process 
onshore-directed  bottom  cu>  nts  may  have  been  generated  giving  rise  to 
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igure  57.  Bottom  current  speed  and  d meet  ion 


Observed  versus  predicted  tidal  cycle  in  the  entrance  channel  for  29  March  through 
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the  NW-trending  vector  path  in  Figure  52  for  the  period  30  March  to  3 
Apri I . 


109.  Progressive  vector  diagrams  are  useful  in  assessing  both 
the  net  flow  direction  and  (assuming  a uniform  flow  field)  the  net  drift 
for  a particular  time  period.  The  following  tabulation  lists  the  calcu- 
lated net  flow  drift  and  direction  for  each  of  the  vector  diagrams 
i I lustrated: 


Location 

Date.  1976 

Duration  of 
Observation 
Period  (Days) 

Net 

Drift  (km) 

Net 

Di recti on 

Buoy  B 

17-24  March 

7 

195 

SW 

8 Apri 1 - 1 May 

24 

417 

sw 

Buoy  D 

9 Oct  - 14  Nov 

37 

260 

SW 

1 7 March  - 7 Apri 1 

22 

390 

SE 

8 April  - 3 May 

26 

277 

SE 

110.  Thus,  for  the  time  periods  of  observation,  the  net  flow 
direction  in  the  study  area  is  oriented  in  an  offshore  or  down-coast 
direction.  Calculated  rates  of  net  drift  range  from  7 to  28  km/day.  The 
implication  of  these  results  to  bedload  transport  of  disposed  dredged 
material  is  apparent.  Assuming  that  the  direction  and  rate  of  net  drift 
are  valid  indicators  of  net  bedload  transport,  then  dredged  material 
located  In  the  disposal  site  will  likely  be  transported  down  coast  or 
offshore  and  thus  not  return  to  the  channel  to  cause  shoaling. 

Oscillatory  bottom  velocities 

111.  Application  of  wave  theory  to  near-bottom  oscillatory 
flow  Indicates  that  such  motion  may  play  a significant  role  In  sediment 
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transport  processes  within  the  offshore  dredged  material  disposal  site. 
Surface  waves  produce  an  oscillatory  motion  which,  near  the  bottom,  reduce 
to  a back-and-forth  horizontal  motion,  the  orbital  diameter  (dQ)  such  that 

dc  = H/sIn  h (kh)  ( 16) 

where  (H)  Is  the  wave  height,  (h)  is  the  water  depth  and  (k)  is  the  wave 
number  that  is  equal  to  2 divided  by  the  wavelength  (Lw).  Wavelength  is 
obtained  by: 

Lw  = Lo  tan  h (2  h/Lp)  (17) 

where  (Lq),  the  deepwater  wavelength,  is  equal  to  I . 56  mu  1 1 i p I led  by  the 

square  of  the  period  in  seconds.  The  maximum  oscillatory  bottom  velocity 
(Um)  attained  by  this  horizontal  bottom  motion  is  then: 

Um  = VT  <lp) 

where  (T)  is  the  wave  period.  The  oscillatory  bottom  velocity  can  be 
predicted  as  a function  of  the  depth  of  water  and  the  wave  height  and 
period. 


112.  Unfortunately,  wave  height  and  period  measurements  in 
the  study  area  were  not  logged  during  field  operations.  However,  annual 
sea  surface  data  were  acquired  from  Tables  18-19  of  "Summary  of  Synoptic 
Meteorological  Observations,"  U.S.  Naval  Weather  Command  (1976).  This 
Information  was  compiled  from  marine  surface  observations  made  by  ships 
In  passage  near  Galveston.  Tables  4-6  are  based  on  a total  of  12,139 
observations  made  during  the  period  |9b3-|97|.  Table  4 shows  the  percen- 
tage of  frequency  of  surface  wind  speed  and  direction  versus  sea  height. 
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Table  4 


Table  4 (continued) 


Wind 

Speed,  m/sec  - 

t 

Heiqht,  m 

0.5-1. 5 

1. 5-5.0 

5.0-10.8 

10.8-17.0 

Total 

Northeast 

3.25-3,75 

.0 

.0 

.1 

.1 

3.75-4,75 

.0 

.0 

.0 

.1 

" : +a  1 

v.ent 

.2 

4.8 

6.8 

1 .4 

13.3 

East 

.25 

.2 

l.l 

.0 

1.3 

.25-. 75 

.2 

3.9 

1.7 

.0 

5.8 

.75-1.25 

2.4 

4.0 

.2 

6.5 

1 .25-1.75 

.5 

2.4 

.3 

3.2 

1 .75-2.25 

.0 

.1 

.8 

.3 

l.l 

2.25-2.75 

.0 

.2 

.2 

.4 

2.75-3.25 

.0 

.0 

. 1 

.1 

3.25-3.75 

.0 

.0 

. 1 

.1 

3.75-4.75 

.0 

.0 

.0 

Total 

F’ercent 

.4 

7.9 

9.1 

1 .0 

18.5 

Southeast 

.25 

.3 

1.5 

.0 

1.8 

.25-. 75 

.1 

4.8 

2.4 

.0 

7.4 

.75-1.25 

.1 

3.0 

5.7 

.2 

9.0 

1.25-1.75 

.6 

3.4 

.3 

4.4 
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Wind  Speed,  m/sec  - 

t 

Heiqht,  m 

0.5-1. 5 

1. 5-5.0 

5.0-10.8 

10.8-17.0 

Total 

Southeast 

1 .75-2.25 

.1 

1 .1 

.3 

4.4 

2.25-2.75 

.3 

.1 

.3 

2.75-3.25 

.0 

.0 

.1 

.1 

.2 

3.25-3.75 

.0 

.0 

3.75-4.75 

.0 

.0 

Tota  1 
Percent 

.9 

10. 1 

13.1 

l.l 

24.8 

South 

.25 

.4 

1 .3 

.1 

.0 

1 .7 

.25-. 75 

.2 

3.6 

1 .6 

.0 

5.3 

.75-1.25 

1 .9 

3.5 

. 1 

5.5 

1 .25-1 .75 

.0 

.3 

1 .9 

.3 

2.5 

1 .75-2.25 

.0 

.1 

.6 

.2 

.9 

2.25-2.75 

.0 

.1 

.2 

.4 

2.75-3.25 

.0 

. 1 

.1 

3.25-3.75 

.0 

.0 

3.75-4.75 

.0 

.0 

.0 

Total 

Percent 

.5 

7.2 

7.8 

.9 

16.9 

Southwest 

.25 

.1 

.5 

.0 

.7 
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Table  4 (continued) 


Wind 

Speed,  m/sec  - 1 

l 

Heiqht.  m 

0.5-1. 5 

1. 5-5.0 

5.0-10.8 

10.8-17.0 

Tota  1 

Southwest 

.25-. 75 

.1 

1.4 

.9 

.0 

1.9 

.75-1.25 

.6 

.7 

1.6 

1.25-1.75 

.0 

.1 

.4 

.5 

1. 75-2. 75 

.0 

.1 

.1 

2.25-2.75 

.0 

.1 

.1 

2.75-3.25 

.0 

.0 

.0 

3.25-3.75 

.0 

.0 

.0 

3.75-4.75 

.0 

.0 

.0 

Total 

Percent 

.2 

2.6 

1.7 

.2 

4 <7 

West 

.25 

.1 

.5 

.0 

.6 

.25-. 75 

.1 

.9 

.3 

.0 

1.3 

.75-1 .25 

.0 

.3 

.4 

.7 

1.25-1.75 

.1 

.3 

.4 

1.75-2.25 

.0 

.1 

.1 

.2 

2.25-2.75 

.0 

2.75-3.25 

.0 

.0 

3.25-3.75 

.0 

.0 

.0 

3.75-4.75 

.0 

.0 

.0 

.0 

.0 

Total 

Percent 

.2 

1.7 

1.2 

.2 

3.3 
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Wind  Speed,  m/sec  - % 


Height,  m 0.5-1 .5  I .5-5.0  5.0-10.8  10.8-17.0 

Northwest 


.25 

.1 

.4 

.0 

•5 

.25-. 75 

.1 

.9 

.3 

.0 

1,2 

.75-1.25 

.4 

.7 

.1 

1.2 

1 .25-1 .75 

.1 

.5 

.2 

.9 

1 .75-2.25 

.0 

.3 

.3 

.7 

2.25-2.75 

.0 

.0 

.2 

.2 

2.75-3.25 

.0 

.0 

.1 

. 1 

3.25-3.75 

.0 

.0 

. 1 

3.75-4.75 

.0 

.0 

. 1 

. 1 

Total 

Percent 

.2 

1.8 

2.0 

1 .0 

VJ1 

O 

TOTAL  PERCENT  =97.6 
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Table  5 


1 

Wind 

Speed  Versus  Sea  Height 

Wind  Speed, 

m/sec 

Height.  m 

0.5-1 .5 

1. 5-5.0 

5.0-10.8 

10.8-17.0 

17.0-24.2 

Total 

.25 

4.1 

6.6 

0.3 

0 

0 

10.9 

.25-. 75 

1.0 
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0 

0 

28.7 

.75-1.25 

0.2 
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0.9 

0 

31.9 

1.25-1.75 

2.0 
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2.1 

0.1 

16.8 

1 .75-2.25 

0.3 

4.3 

2.2 

0.1 

6.9 

2.25-2.75 

0.1 

1.2 

1.4 

0.1 

2.9 

2.75-3.25 

0 

0.3 

0.7 

0.1 

1 .1 

3.25-3.75 

0 

0.1 

0.3 

0.5 

3.75-4.75 

0 

0 

0.1 

0.2 

0 

0.3 

Total 

Percent 

5.3 

39.1 

47.0 

7.8 

0.4 

100.0 

L , A 


Percent  Frequency  of  Wave  Height 


Table  5 contains  wind  speed  versus  sea  height.  Table  6 lists  the  percent 
frequency  of  wave  height  versus  wave  period.  Important  findings  gained 
from  an  examination  of  these  tables  are: 


a.  Prevailing  wind  azimuths  for  offshore  Galveston  in 
decreasing  order  of  frequency  are  southeast,  east,  and 
south. 

b.  Common  wind  speeds  ranging  from  1.5  to  10.8  m/sec 
generate  sea  waves  0.25  to  1.75  m in  height. 

c.  Surface  waves  ranging  in  height  from  0.25  to  1.75  m 
usually  have  wave  periods  less  than  6.0  sec.  However, 
swells  exceeding  1.25  m in  amplitude  have  occasion- 
ally been  recorded  with  periods  ranging  from  6.0  to 
9.0  sec  in  duration. 

113.  The  significance  of  these  results  in  assessing  the  effect 
of  oscillatory  motion  on  bottom  sediment  is  revealed  by  a plot  of  the 
oscillatory  bottom  velocity  versus  the  wave  period  for  an  average  14.0-m 
water  depth  (Figure  60).  This  graph  was  plotted  using  Equations  16,  17, 
and  18.  As  shown  here,  wave  periods  in  excess  of  5.0  sec  in  duration 
coupled  with  wave  heights  greater  than  0.75  m theoretically  induce  oscil- 
latory bottom  speeds  exceeding  10  cm/sec.  The  percent  frequency  of  wave 
height  versus  wave  period  data  in  Table  6 suggest  that  such  wave  ampli- 
tudes and  periods  can  occur  up  to  30  percent  of  the  time.  If  these 
oscillatory  bottom  velocities  are  superimposed  upon  an  existent  10-  to 
20-cm/sec  bottom  current,  then  the  threshold  of  silt  and  clay  movement 
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( t ^ * 1.0  dyne  /cmM  will  probably  be  exceeded.  Such  movement  may 
then  allow  particles  to  be  entrained  within  a given  flow.  Continuous 
near-bottom  oscillatory  motion  may  also  prevent  deposition  of  suspended 
particles  If  the  resultant  bed  shear  exceeds  the  threshold  shear  stress 
dyne  /cm*)  permitting  deposition.  Hence,  sea-state  conditions  in 
the  dredged  material  disposal  site  can  have  a significant  effect  upon 
local  near-bottom  sediment  transport  processes. 

Summary 

114.  The  physical  oceanograph i c and  meteoro logic  data  acquired 
during  the  period  February  1^75  through  May  1*776  for  offshore  Galveston 
and  vicinity  indicate  that: 

a.  Tides  in  the  entrance  channel  are  mixed  in  nature  and 
typically  range  from  0.3  to  0.6  m in  height. 

b.  Tidal  currents  originating  from  cyclic  passage  through 
the  entrance  channel  induce  reversing  rotary  circula- 
tion patterns  over  the  offshore  dredged  material 
disposal  site.  In  particular,  ebb  tide  flows  generate 
S-SW  directed  currents,  whereas  flood  tide  currents 
produce  E-NE  directed  flows  during  maximum  current 
development. 

c . A shear  zone  periodically  exists  near  the  bottom  of  the 
water  column  in  the  offshore  disposal  are.  This  layer 
appears  to  form  in  response  to  either  reversing  tidal 
currents  or  near-bottom  circulation  patterns  set  up  by 
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prevailing  storm-wind  azimuths. 


d . Adverse  meteorological  conditions  affect  both  the 
nearshore  sea  level  datum  and  the  near-bottom  circula- 
tion pattern.  In  particular,  primary  storm-wind 
azimuths  affecting  the  hydrographic  regime  include: 

(I)  strong  persistent  E,  SE,  or  S winds  traversing 
the  Gulf  of  Mexico  and  (2)  strong  N to  NW  winds 
following  the  passage  of  a cold  front 

e.  Maximum  bottom  current  velocities  and  net  offshore- 
directed  to  down-coast  bottom  flows  appear  to  be 
generated  by  strong  prevailing  southeasterly  winds 
in  excess  of  about  5 m/sec. 

f . Theoretical  consideration  of  near-bottom  oscillatory 
flows  induced  by  surface  waves  suggests  that  such 
motion  may  play  a significant  role  in  governing  near- 
bottom sediment  transport  processes,  particularly 
during  adverse  meteorological  conditions. 

Applicability  of  Flume  Experimental  Results 
to  Offshore  Sediment  Transport  at  Galveston 

Scaling  considerations 

115.  The  applicability  of  the  results  obtained  from  these  flume 
experiments  to  offshore  sediment  transport  can  be  partially  realized  if 
consideration  is  given  to  scaling  effects  between  the  prototype  and  the 
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natural  environment.  Dynamic  similitude  between  the  open  channel  flow 
and  the  offshore  currents  requires  reproduction  of  primarily  the  gravity, 
inertial,  and  viscous  forces  that  characterize  these  ocean  currents. 
Unfortunately,  there  does  not  exist  a single  relationship  that  combines, 
as  ratios,  all  three  of  these  forces.  However,  scaling  on  the  basis  of 
the  Reynolds  number  for  both  full  and  scaled  conditions  preserves  a 
constant  ratio  between  the  viscous  and  inertial  forces.  Similarly, 

Froude  number  (Fr)  scaling  is  based  upon  a constant  ratio  between  gravity 
and  inertial  forces.  Thus,  these  two  scaling  parameters  were  used  to 
ascertain  those  open  channel  flow  conditions  that  displayed  dynamic 
similarity  to  the  ocean  environment. 

116.  Appropriate  values  of  the  Reynolds  number,  as  previously 
discussed  (refer  to  paragraph  16),  establish  whether  a given  flow  is  tur- 
bulent or  laminar  in  nature.  Due  to  the  inherently  large  length  scale 
needed  to  characterize  ocean  currents,  the  Reynolds  number  will  result 
in  values  exceeding  several  thousand  for  current  speeds  greater  than  a 
few  centimeters  per  second  (see  Equation  5).  Such  currents  are  typically 
considered  as  fully  turbulent  in  nature.  Reynolds  number  calculation  for 
open  channel  flow  in  the  flume  (see  Appendix  E)  indicated  that  turbulent 
flow  (R0  > 2000)  theoretically  occurred  at  current  velocities  in  excess 
of  about  0.3  cm/sec.  In  addition,  observation  of  fluorescent  dye  injected 
into  the  observation  channel  water  column  during  flow  rates  ranging  from 
5 to  10  cm/sec  showed  full  turbulent  eddy  development  with  time  downstream 
from  the  injection  point.  Thus,  the  existence  of  fully  developed  turbulent 
flow  during  experimental  runs  was  confirmed. 
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117.  the  Froude  number  describes  the  relationship  between  the 


Inertial  and  gravity  forces  for  a particular  flow  such  that: 

U 


where  (L)  represents  the  reference  lenoth  and  (of  Is  the  accelernt Ion  due 
to  gravity.  Calculation  of  egulvalont  Troude  and  Reynolds  numbers  between 
the  prototype  and  observed  bottom  currents  (Appendix  I ) established  that 
dynamic  similarity  in  the  vertical  existed  In  the  open  channel  flow  for 
current  speeds  up  to  about  cm/sec.  This  similarity  exists  for  uniform 
flow  conditions  that  are  characterized  by  fully  developed  turbulent  eddies. 

Critical  erosion  and  deposition  velocities 

lift.  Application  of  the  experimentally  determined  critical 
erosion  and  deposition  velocities  as  measures  of  the  critical  velocities 
occurring  within  the  offshore  disposal  slfe  Is  valid  under  certain  assump- 
tions. These  assumptions  Include  the  following:  (I)  bottom  currents  In 
the  offshore  dredged  material  site  are  characterized  by  uniform,  fully 
developed  turbulent  flow;  (,?)  the  effect  of  wave-induced  oscillatory 
motion  on  the  threshold  of  cohesive  sediment  movement  is  negligibly  small; 
and  (')  the  hydrodynamlcal  roughness  exhibited  by  the  dredged  material 
placed  into  the  flume  is  similar  to  the  roughness  existent  on  the  offshore 
disposal  mounds,  tstab I I shment  of  the  validity  and/or  variability  of  those 
three  assumptions  requires  the  development  of  new  field  monitoring  technique' 
to  measure  the  in  situ  bottom  current  and  sediment  conditions  over  long  time 
periods.  Until  such  techniques  are  developed,  extrapolation  of  the  "criti- 
cal velocity"  results  to  field  conditions  should  he  viewed  as  a best 
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estimate  based  upon  state-of-the-art  methods. 


Sediment  Transport  Processes  in  the  Offshore 
Dredged  Material  Disposal  Site 


E t feet  of  c I imate 

{19.  Persistent  winds.  Examination  of  both  metoorologlc  and 
hydrographic  data  tor  offshore  Galveston  during  the  period  February  |97*s 
through  May  |P7(>  indicates  the  importance  of  local  climate  on  sedimento- 
loglc  processes  operative  in  the  marine  environment.  In  particular,  the 
intensity  and  duration  of  persistent  winds  Influence  the  nature  and  direc- 
tion ot  offshore  bottom  currents  that  erode,  transport,  and  deposit  bottom 
sediments. 


120.  The  principal  wind  azimuths  affecting  the  Galveston 
vicinity  are  illustrated  as  a surface  wind  rose  in  Figure  61.  Those  wind 
directions  represent  the  relative  percent  frequency  of  prevailing  wind 
directions  as  listed  in  Table  4 for  the  period  1 965-1971,  Also  shown  In 
Figure  61  are  net  bottom  drift  vectors.  The  longshore-d I rected  vectors 
located  adjacent  to  the  Galveston  Island  shoreline  and  those  within  the 
entrance  channel  are  based  upon  the  sediment  transport  data  presented  by 
Moll  (1976).  The  net  bottom  drift  vectors  Illustrated  at  Puoy  IF  and  Puoy 
D are  positioned  according  to  the  net  drift  directions  given  by  the 
progressive  vector  trajectories  illustrated  In  Figures  60  - 64. 

121.  The  net  drift  azimuths  shown  at  Puoys  P and  P presumably 
resulted  from  strong  prevailing  southeasterly  winds.  It  should  be 
realized  that  these  bottom  drift  indicators  are  representative  only  of  the 
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periods  during  which  the  bottom  current  meters  were  operative.  This  obser- 
vation period,  which  totaled  31  days  at  Buoy  B and  85  days  at  Buoy  D,  is 
not  extensive  enough  in  duration  to  establish  annual  net  drift  directions 
with  statistical  confidence.  However,  based  upon  the  correlation  found 
between  the  persistent  southeasterly  winds  and  net  longshore  to  offshore 
bottom  flow  azimuths,  it  appears  reasonable  to  speculate  that  annual  net 
drift  directions  are  oriented  similarly.  This  speculation  is  based  upon 
the  fact  that  southeasterly  winds  occur  most  frequently  on  a yearly  basis 
as  evidenced  by  the  8-yr  averaged  prevailing  wind  data  in  Table  4. 

122.  Tropical  storms.  A discussion  of  climatic  effects  upon 
sediment  transport  processes  for  offshore  Galveston  would  not  be  complete 
without  consideration  of  severe  tropical  storms.  Although  no  hurricanes 
seriously  affected  the  Galveston  vicinity  during  the  period  of  this  study, 
destructive  storms  such  as  hurricane  Audrey  (1957)  and  Carla  (1961)  have 
produced  extreme  storm  tides.  During  the  passage  of  Audrey,  a positive 
tide  anomaly  of  3.1  m was  recorded  at  Galveston.  Similarly,  Carla  caused 
a 2.3-m  tidal  height  rise.  Hayes  and  Scott  (1964)  documented  the  effects 
of  the  storm  surge  and  associated  currents  generated  by  hurricane  Carla 
along  the  central  Texas  coast.  After  passage  of  the  storm  inland,  seaward- 
flowing density  currents  deposited  a fresh  sand  layer  (up  to  3 cm  thick) 
over  a sandy  mud  bottom  in  depths  of  15  to  20  m.  In  addition,  a graded 
bed  (up  to  10  cm  thick)  comprised  of  very  fine  sand,  silt,  and  clay  was 
traced  over  the  homogenous  mud  bottom  to  depths  exceeding  60  m.  Such 
sedimento logic  features  point  out  the  major  role  that  a tropical  storm  can 
have  in  modifying  an  offshore  geologic  environment  over  a very  short  time 
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Modes  of  sediment  transport 

123.  Information  on  the  various  modes  of  sediment  transport 
for  offshore  Galveston  was  obtained  by  combining  the  results  of  the 
flume  experiments  with  those  obtained  from  the  hydrographic  data  in  this 
study.  The  mo^t  valuable  field  data  used  to  assess  the  relative  degrees 
of  erosion,  transport,  and  deposition  occurring  in  the  disposal  site  were 
the  current  speed  and  direction  measurements  recovered  from  current  meters 
positioned  at  Buoy  B and  Buoy  D I m above  the  bottom.  Discrete  15-min 
interval  current  vector  readings  were  analyzed  for  the  relative  frequency 
of  occurrence  of  cumulative  5.0-cm/sec  current  speed  increments  over  the 
total  range  of  values  recorded  at  each  station.  The  flow  velocity  distri- 
bution of  each  observation  period  for  both  stations  was  then  determined. 
Cumulative  curves  showing  the  percentage  distribution  of  flow  speed  incre- 
ments (plotted  on  a probability  axis)  versus  the  logarithm  to  the  base  ten 
of  the  measured  current  speeds  (plotted  on  an  arithmetic  axis)  served  as 
the  basis  for  these  determinations.  These  curves  are  illustrated  in  Figures 
62  through  66.  The  velocity  distribution  shown  in  Figure  62  was  plotted 
from  cumulative  speed  data  given  by  Hall  (1976). 

124.  Examination  of  these  velocity  distributions  indicates 
that  the  median  velocities  displayed  by  the  curves  for  Buoy  B ranged  from 
22  to  35  cm/sec,  while  those  for  Buoy  D were  slightly  lower,  ranging  from 
18  to  22  cm/sec.  However,  another  purpose  of  presenting  the  current  speed 
data  for  a specific  time  period  in  the  form  of  a cumulative  curve  was  to 
ascertain  the  relative  frequency  during  which  erosion,  transportation,  and 
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Figure  62.  Cumulative  curve  of  the  current  velocity  distribution  for 
9 October  to  14  November  1975  at  Buoy  D modified  from 
Hall  (1976)  (0  = average  current  speed  measured  I m 
above  bed) 
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Figure  63.  Cumulative  curve  of  the  current  velocity  distribution 
for  17  to  24  March  1976  at  Buoy  D (0  = average  current 
speed  measured  I m above  bed) 
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Figure  64.  Cumulative  curve  of  the  current  velocity  distribution 
for  17  March  to  7 April  1976  at  Buoy  D (0  = average 
current  speed  measured  I m above  bed) 
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Figure  66.  Cumulative  curve  of  the  current  velocity  distribution 

for  8 April  to  3 May  1976  at  Buoy  0 (0  = average  current 
speed  measured  I m above  bed) 


deposition  occurred  at  a specific  locality.  The  frequency  of  the  current 
speeds  for  each  observation  period  was  examined  in  order  to  determine  the 
percent  of  the  time  that  the  critical  erosion  velocity  was  exceeded.  This 
percentage  then  provided  a relative  assessment  of  the  time  that  rapid  bed- 
load erosion  may  have  occurred  at  a particular  locality.  Likewise,  the 
frequency  at  which  current  speeds  occurred  below  the  critical  deposition 
velocity  Indicated  a relative  estimate  of  the  time  during  which  perhaps 
significant  tine  silt  and  clay  deposition  occurred.  Finally,  the  fre- 
quency of  speeds  existent  between  the  critical  erosion  and  deposition 
velocities  approximated  the  time  during  which  mostly  transport  of  sus- 
pended sediment  likely  occurred. 

125.  The  following  tabulation  lists  the  relative  percentage 


frequency  of  erosion,  transport,  and  deposition  that  occurred  at  Buoy  B 
and  Buoy  D per  given  observation  period: 


Duration  of 

DEPOSITION 

TRANSPORT 

EROSION 

Date  of 

Observation 

% of  Time 

% of  Time 

i of  Time 

Observation 

Period,  days 

U < 1 2 cm/sec 

Buoy  B 

1 2 < U < 28 

0 > 28  cm/ sec 

3/17/76- 

3/24/76 

7 

3 

16 

81 

4/  8/76- 
5/  1/76 

24 

3 

Buoy  D 

44 

53 

10/  9/75- 
11/14/75 

37 

19 

59 

22 

3/17/76- 
4/  7/76 

22 

20 

49 

31 

4/  8/76- 
5/  3/76 

26 

18 

56 

26 

The  values  of  the  critical  erosion  and  deposition  velocities  (U  = 28  and 
12  cm/sec,  respectively)  shown  in  the  tabulation  were  ascertained  from 
Einstein's  Bedload  Function  for  open  channel  flow  over  smooth  boundaries 
using  the  experimentally  determined  critical  erosion  and  deposition  bed 

shears  (t  =1.0  and  0.2  dyne  /cm^,  respectively)  to  calculate  values 

oc  r 

of  the  friction  velocity.  Examininq  the  values  in  the  tabulation,  it  is 
apparent  that  rapid  cohesive  bedload  erosion  probably  occurred  much  more 
frequently  at  Buoy  B than  at  0.  Similarly,  current  speeds  permitting 
significant  silt  and  clay  deposition  likely  occurred  more  often  at  Buoy 
D than  at  Buoy  B.  The  relative  percentage  of  time  that  transport  condi- 
tions (12  < 0 < 28  cm/sec)  occurred  at  these  sites  does  not  necessarily 
indicate  the  total  time  during  which  a sedimentary  particle  remained  in 
suspension.  For  instance,  if  a silt  or  clay  particle  was  deposited  from 
suspension  on  a cohesive  bed  during  a decreasing  flow  rate,  such  as  during 
slack  water  development  during  a tide  cycle,  that  particle  would  probably 
adhere  to  the  bed  until  the  bottom  velocity  again  Increased  to  the  criti- 
cal erosion  velocity  needed  to  resuspend  it.  The  bottom  current  may  not 
resume  erosional  competency  until  maximum  ebb  or  flood  tidal  current 
development  is  reinstated  several  hours  later.  In  addition,  both  erosion 
and  deposition  of  fine  silt  and  clay  can  occur  during  the  range  of  these 
so-called  ''transport”  velocities  as  evidenced  by  ( I ) the  negative  and 
positive  slopes  of  the  scour  curves  for  the  dense  and  floe  bed  erosional 
runs  and  (2)  the  observed  decrease  in  the  silt  content  (increase  in  the 
clay  content)  of  the  washload  during  sustained  flow  rates.  Thus,  these 
critical  erosion  and  deposition  velocities  are  nonunique  and  delineate 

only  the  relative  degree  of  erosion,  transport,  or  deposition  that 
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occurred  at  a given  locality. 


Offshore  disposal  mounds 

126.  The  relative  degrees  of  erosion  and  deposition  occurring 
at  each  of  the  disposal  mounds  (see  Figure  2)  should  be  indicated  indirectly 
by  the  texture  of  the  surflcial  sediments,  and  also  by  general  changes  in 
the  bathymetric  configuration  of  the  mounds.  For  instance,  surficial 
dredged  material  at  Buoy  B would  be  expected  to  display  a coarser  sediment 
texture  than  at  Buoy  D due  to  the  greater  frequency  of  erosional  bottom 
current  velocities  occurring  at  Buoy  B.  Conversely,  the  silt  and  clay 
content  of  surficial  sediment  at  Buoy  D should  be  hioher  than  Buoy  B due 
to  the  greater  frequency  of  depositional  current  velocities  occurring  at 
Buoy  D. 


127.  Expected  variations  in  the  bathymetric  configuration  over 
time  Include  a lower  topographic  relief  at  Buoy  B than  at  Buoy  D (or  Buoy 
C)  and  general  changes  in  the  mound  orientation.  The  lower  topographic 
relief  expected  at  the  Buoy  B mound  is  again  predicated  on  the  evidence 
that  a higher  frequency  of  erosive  bottom  currents  occurs  at  Buoy  B than 
at  D.  Thus,  local  bedload  erosion  would  result  in  the  removal  of  more 
material  at  B than  at  C or  D. 

128.  Finally,  a southwest-  and/or  southeast-trending  mound 
configuration  might  be  expected  to  form  at  these  disposal  mounds  if  the 
net  bottom  drift  directions  shown  in  Figure  61  are  valid  indicators  of 
net  bed  load  movement.  The  analysis  of  variations  in  the  bottom  sediment 
textural  distributions  and  bathymetric  mound  configurations  as  a function 
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of  time  by  Estes  and  Scruda+o  (1977)  may  indicate  such  general  trends. 


PART  V:  CONCLUSIONS 


Flume  Experiments 

129.  As  a result  of  the  flume  experiments  performed  on  the 
four  sand,  silt,  and  clay  mixtures,  it  was  concluded  that: 

a.  Each  of  the  four  sediments  eroded  similarly  as  evi- 
denced by  the  similar  critical  bed  shear  value 

(x  r **  I .0  dyne  /cm^)  necessary  to  initiate  massive 
bedload  transport  and  rapid  TSM  concentration 
I ncroases. 

b.  Rapid  deposition  of  suspended  fine  silt  and  clay  occurred 

2 

at  very  low  bed  shear  values  (t  < 0.2  dyne  /cm  ). 

c.  Primary  silt  and  clay  resusponsion  occurred  durlna  the 
initial  2 to  4 hr  fol lowing  each  current  velocity  In- 
crease. Thereafter,  for  sustained  flow  rates,  the 
scour  rate  decreased  to  a near  constant  level . 

d.  Continuous  erosion  on  a hydrodynamlcal ly  rough  bed 
resulted  in  textural ly  diverse  sediment  interfaces 
within  the  sediment  column. 

e.  Washload  composition  consisted  primarily  of  fine  silt 
and  clay  at  a 30-cm  elevation  above  the  bed.  Signifi- 
cant variations  in  the  washtoad  grain-size  distribu- 
tion occurred  with  Increaslnq  or  decreasing  flow  rates. 
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Hydrography  of  Pffshoro  Galveston 

130.  It  is  concluded  from  an  analysis  of  the  hydrographic 
data  obtained  for  offshore  Galveston  that: 

a.  Tidal  currents  produced  by  cyclic  passage  through 
the  Galveston  Bay  entrance  channel  Induce  reversing 
rotary  circulation  patterns  over  the  offshore  dredged 
material  disposal  site.  In  particular,  ebb  tide  flows 
generate  S-SW  directed  currents,  whereas  flood  tide 
flows  generate  E-NE  directed  currents. 

b . A shear  zone  periodically  exists  near  the  bottom  of  the 
water  column  in  the  offshore  disposal  area.  This  zone 
appears  to  form  in  response  to  either  reversing  tidal 
currents  or  near-bottom  circulation  patterns  set  up  by 
prevailing  storm  wind  azimuths. 

c . Adverse  meteorological  conditions  affect  both  the 
coastal  sea  level  datum  and  the  near-bottom  circu- 
lation pattern.  Primary  storm  wind  azimuths  affecting 
the  hydrographic  regime  Include:  (I)  strong  per- 
sistent southeast,  east,  or  south  winds  traversing 
across  the  Gulf  of  Mexico  and  (2)  strong  north  to 
northwest  winds  following  the  passage  of  a cold  front. 

d . Maximum  bottom  current  velocities  and  net  down-coast 
to  offshore-directed  bottom  flows  appear  to  be 
generated  by  southeaster  I y storm  winds  exceeding 

Ibl 

L 


about  b m/sec. 


Offshore  Sediment  Transport 

131.  Based  on  the  evaluation  of  both  the  offshore  hydrographic 
and  laboratory  experimental  results.  It  appears  that: 

a . Erosion  Is  occurring  more  frequently  near  the 
northern  margin  (near  Buoy  B)  of  the  offshore  dis- 
posal site  than  at  the  southern  margin  (near  Buoy 
D).  Conversely,  deposition  of  suspended  silt  and 
clay  occurred  more  frequently  near  the  southern 
boundary  than  at  the  northern  end. 

b.  Net  bed  load  transport  of  dredged  material  moves 

down  coast  or  offshore  from  the  Galveston  Bay  entrance 
channel.  Thus,  the  probability  of  significant  quan- 
tities of  dredged  material  returning  to  the  entrance 
channel  as  shoaling  sediment  is  small. 
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APPENDIX  A:  TOTAL  SUSPENDED  MATTER  CONCENTRATION,  TIME 
AFTER  VELOCITY  CHANGE,  MEAN  FLOW  SPEED,  AND  WATER 
LOCATIONS  OBTAINED  DURING  EXPERIMENTAL  RUNS 


Table  A1 


Total  Suspended  Matter  Concentration.  Time  After 


Velocity  Chanqe.  Mean  Flow  Speed,  and  Water  Sample 


Locations  Obtained  Durlna  Experimental  Runs 


Set  Run 


Sample  # 


Ut 

(cm/sec) 


Tvtt  TSMt 


(mg/  f ) 


Buoy  C 1 

1 

02 

30.5 

10 

0.1 

0.95 

2 

02 

0.5 

1 .42 

3 

E2 

0.5 

1.64 

« 

4 

02 

1 .0 

2.07 

5 

E2 

1 .0 

3.00 

6 

C2 

2.0 

4.97 

1 

7 

02 

2.0 

5.28 

j 

8 

°2 

4.0 

7.12 

12 

°2 

16.0 

7.97 

14 

02 

24.0 

7.86 

Buoy  C 2 

15 

d2 

30.5 

15 

0. 1 

8.17 

f 

17 

D2 

1.0 

9.56 

18 

o2 

2.0 

12.25 

20 

d2 

6.0 

14.38 

22 

°2 

12.0 

17.51 

23 

02 

22.3 

18.67 

25 

d2 

24.0 

17.60 

Buoy  C 3 

26 

d2 

30.5 

25 

O'.  1 

17.98 

27 

0.5 

24.55 

28 

0.9 

26.13 

32 

8.0 

31 .05 

*LS  = Sample 

location. 

see 

Figure 

19. 

**Z  = Water  depth  of  sample. 

fU  = Mean  flow  speed. 

t!Tv  = Time 

after  velocity 

change. 

tTSM  = Total 

suspended 

matter  concentration. 

' 

A2 
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Table  A1  (continued) 


Set  Run 

Sample  If 

Ls 

Z 

(cm) 

U 

(cm/sec) 

Tv 

(hr) 

TSM 
(mg// ) 

Buoy  C 3 

34 

d2 

30.5 

23 

22.3 

37.34 

Buoy  C 4 

37 

d2 

30.5 

31 

0.1 

29.14 

38 

0.5 

59.75 

39 

1 .0 

67.00 

40 

2.0 

80.19 

41 

4.0 

83.75 

42 

6.2 

95.36 

43 

8.0 

90.59 

44 

12.0 

95.99 

45 

13.3 

90.71 

48 

35.4 

95.05 

49 

47.4 

89.53 

50 

59.0 

96.75 

51 

71.6 

94.82 

52 

82.3 

99.49 

53 

98.4 

98.73 

Buoy  C 5 

54 

02 

30.5 

23 

0.1 

91.54 

55 

E2 

0.2 

94.14 

56 

D2 

0.3 

98.08 

57 

°2 

7.0 

93.09 

58 

02 

8.4 

84.18 

60 

o2 

23.4 

88.48 

Buoy  C 6 

61 

d2 

30.5 

10 

0.1 

89.56 

62 

0.7 

93.18 

63 

8.5 

81.63 

64 

10.3 

80.00 

65 

19.9 

78.45 

Block  15  2 

74 

d2 

30.5 

16 

0.1 

3.06 

75 

1.9 

2.52 

77 

8.3 

4.26 

79 

15.3 

5.58 

80 

25.5 

4.76 

Block  15  3 

81 

02 

30.5 

23 

0.1 

9.46 

82 

0.7 

6.62 

83 

7.0 

19.57 

87 

1 1.0 

21 .85 

89 

24.8 

19.62 

A3 


Table  A1  (continued) 


Set  Run 

Sample  0 

Ls 

Z 

(cm) 

0 

(cm/sec) 

Tv 
( hr) 

TSM 
(mg/f ) 

Block  14  4 

90 

d2 

30.5 

31 

. 1 

73.54 

91 

2.75 

150.12 

92 

6.25 

182.52 

95 

10.0 

163.67 

97 

23.0 

175.66 

Block  15  5 

98 

°2 

30.5 

10 

.1 

75.52 

101 

6.8 

72.02 

102 

1 1 .0 

88.18 

103 

22.8 

63.38 

Block  27  1 

105 

C2 

30.5 

12 

.3 

15.63 

106 

E2 

1 .7 

20.28 

107 

E2 

2.8 

19.83 

108 

e2 

8.4 

20.79 

109 

e2 

46.5 

15.88 

1 10 

E2 

51.0 

15.65 

Block  27  2 

1 1 1 

E2 

30.5 

23 

. 1 

14.47 

112 

.2 

19.63 

113 

.5 

24.49 

1 14 

1 .0 

29.91 

1 15 

1 .5 

1 16 

2.0 

34.88 

1 17 

4.0 

36.45 

1 18 

6.0 

36.66 

1 19 

9.1 

35.74 

120 

11.9 

3*.8I 

121 

24.9 

33.57 

Block  27  3 

122 

e2 

30.5 

30 

.2 

36.57 

123 

.5 

37.50 

124 

1.0 

43.58 

125 

1.5 

46.53 

126 

2.0 

48.57 

127 

4.0 

49.28 

128 

6.0 

45.90 

129 

9.3 

45.66 

130 

14.8 

52.17 

131 

22.9 

59.78 

A4 


Table  A1  (continued) 


Set  Run  Sample  # 


Block  27  4 132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

Block  27  5 145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

Block  27  6 157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 


2 U 

L3  (cm)  (cm/sec) 


E2  30.5  38 


E2  30.5  44 


E2  30.5  57 


E 


2 


( hr) 

TSM 
(mg/f ) 

.1 

56.94 

.4 

57.20 

1 .4 

66.91 

1.8 

65.45 

2.0 

73.97 

4.0 

81 .07 

6.6 

89.37 

9.0 

90.13 

12.0 

92.39 

23.9 

102.27 

28.2 

106.04 

37.8 

108.00 

48.0 

109.58 

.1 

115.43 

.5 

1 17.99 

1 .0 

1 16.65 

1.8 

128.00 

4.0 

128.25 

5.5 

149. 1 1 

5.5 

146.18 

5.5 

142.34 

9.8 

155.70 

12.1 

158.04 

23.6 

179. 10 

24.0 

179.88 

.1 

173.70 

.5 

207.89 

1.3 

237.82 

1 .8 

248.56 

6.4 

398.66 

9.0 

435.62 

13. 1 

471.00 

25.6 

670.68 

1 .6 

646.23 

3.8 

643.43 

9.7 

645.74 

23.6 

614.30 

Block  27  7 


30.5 


30 


Table  A1  (continued) 


Set 

Run 

Sample  If 

Ls 

Z 

(cm) 

0 

(cm/sec) 

Tv 

(hr) 

TSM 
(mg/ l) 

Buoy  D 

2 

69 

d2 

30.5 

23 

.1 

29.89 

70 

d2 

30.5 

.5 

36.09 

71 

d2 

15.2 

1 .0 

39.24 

72 

d2 

30.5 

1.0 

40.35 

73 

D2 

45.7 

1.0 

40.00 

74 

d2 

30.5 

1.5 

42.75 

75 

°2 

30.5 

2.0 

40.32 

76 

°2 

30.5 

3.9 

45.60 

77 

d2 

15.2 

6.0 

68.11 

78 

d2 

30.5 

6.0 

45.97 

79 

D2 

45.7 

6.0 

47.00 

80 

d2 

30.5 

8.8 

46.81 

81 

d2 

30.5 

11.7 

44.69 

82 

d2 

15.2 

23.0 

41.36 

83 

d2 

30.5 

23.0 

43.43 

87 

d2 

45.7 

23.0 

43.37 

Buoy  0 

3 

88 

d2 

30.5 

32 

. 1 

42.93 

89 

d2 

30.5 

.5 

61.43 

98 

d2 

30.5 

1 .3 

71  .47 

99 

d2 

15.2 

1.6 

74.36 

100 

d2 

30.5 

1 .6 

74. 1 1 

101 

d2 

45.7 

1 .6 

72.76 

102 

d2 

30.5 

2.0 

69.04 

103 

02 

30.5 

4.1 

70.92 

104 

02 

15.2 

6.0 

71.32 

106 

02 

30.5 

6.0 

69.57 

107 

d2 

45.7 

6.0 

81 .75 

109 

°2 

30.5 

9.1 

80.13 

1 10 

d2 

30.5 

12.0 

77.03 

1 1 1 

°2 

15.2 

23.9 

75.33 

1 12 

d2 

30.5 

23.9 

77.96 

113 

d2 

45.7 

23.9 

77.63 

Buoy  D 

4 

114 

d2 

30.5 

39 

.2 

78.82 

1 16 

D| 

15.2 

.6 

81 .65 

117 

°2 

30.5 

.6 

83.88 

118 

03 

45.7 

.6 

77.87 

1 19 

d2 

30.5 

1 .0 

88.02 

120 

d2 

30.5 

1 .6 

93.78 

121 

d2 

30.5 

2.0 

91.17 

122 

d2 

30.5 

4.0 

99.60 

A7 


Table  A1  (continued) 


Set 

Run 

Sample  # 

Ls 

Z 

(cm) 

0 

(cm/sec) 

(hr) 

TSM 
(mg// ) 

Buoy  D 

4 

123 

Dl 

15.2 

39 

6.3 

100.95 

124 

d2 

30.5 

6.3 

100.24 

12b 

03 

45.7 

6.3 

103.14 

126 

d2 

30.5 

9.0 

103.27 

127 

d2 

30.5 

12.9 

104.79 

219 

0| 

15.2 

23.7 

105.52 

130 

d2 

30.5 

23.7 

102.90 

131 

d3 

45.7 

23.7 

104.93 

Buoy  0 

5 

132 

d2 

30.5 

44 

.1 

105.53 

133 

Dl 

15.2 

.5 

107.42 

134 

d2 

30.5 

.5 

104.19 

135 

03 

45.7 

.5 

|03.46 

136 

02 

30.5 

1.0 

1 10.97 

137 

02 

30.5 

1.5 

107.50 

138 

02 

30.5 

2.0 

106.24 

140 

°2 

30.5 

4.0 

105.11 

141 

°2 

30.5 

6.0 

111.22 

142 

Dl 

15.2 

9.0 

113.35 

143 

02 

30.5 

9.0 

1 14.59 

144 

°3 

45.7 

9.0 

1 16.33 

145 

°2 

30.5 

12.0 

117.17 

146 

Dl 

15.2 

23.0 

119.22 

147 

°2 

30.5 

23.0 

115.55 

148 

03 

45.7 

23.0 

117.55 

Buoy  D 

6 

149 

o2 

30.5 

60 

.1 

125.73 

150 

D| 

15.2 

.5 

164.15 

152 

02 

30.5 

.5 

164.59 

153 

03 

45.7 

.5 

166.74 

154 

02 

30.5 

1.0 

185.14 

155 

o2 

30.5 

1 .5 

192.35 

156 

d2 

30.5 

2.0 

193.60 

157 

02 

30.5 

4.0 

217.70 

158 

D| 

15.2 

6.1 

227.85 

159 

02 

0.5 

6.1 

226.77 

160 

03 

45.7 

6.1 

227.86 

161 

02 

30.5 

9.0 

224.45 

163 

o2 

30.5 

12.0 

212.11 

1 1>4 

D| 

15.2 

24.0 

241.10 

165 

02 

30.5 

24.0 

239.37 

166 

°3 

45.7 

24.0 

223.00 

A8 


Sat  Run  Sample  0 


Buoy  D 


Buoy  0 


Buoy  0 


Buoy  D 1 0 


Buoy  D I I 


(cm/ sec) 


TSM 
(mg// ) 


168 

02 

30.5 

20 

.1 

214.42 

169 

D| 

15.2 

1 .4 

21 1 .72 

170 

d2 

30.5 

1 .4 

213.03 

172 

03 

45.7 

1 .4 

221 .77 

176 

02 

30.5 

5.9 

231 .10 

177 

02 

30.5 

9.4 

230.00 

179 

d2 

30.5 

13.2 

217.41 

180 

°2 

30.5 

23.9 

206.90 

181 

°2 

30.5 

4 

.1 

183.89 

183 

.5 

199.36 

184 

1 . 1 

186.64 

186 

1 .5 

176.32 

1 

4.0 

121.13 

2 

6.1 

117.88 

3 

10.9 

89.87 

4 

13.6 

79.38 

5 

26.0 

53.61 

6 

o2 

30.5 

16 

• *- 

39.31 

7 

.6 

35.61 

0 

1 .0 

34.97 

10 

1.5 

35.44 

1 1 

2.0 

34.97 

187 

4.4 

56.63 

188 

9.0 

29.12 

189 

12.0 

30.92 

190 

24.0 

30.78 

191 

°2 

30.5 

24 

.1 

34.71 

192 

.5 

40.00 

193 

1 .0 

39.58 

194 

1 .6 

42.10 

195 

2.0 

43.62 

196 

4.0 

43.77 

197 

6.2 

46.58 

198 

12.5 

47.  10 

199 

24.5 

48.93 

200 

02 

30.5 

31 

.1 

51  .66 

201 

0? 

30.5 

.5 

59.03 

202 

0| 

15.2 

1 .0 

68.83 
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APPENDIX  B:  CURRENT  VELOCITY  PROFILES  AND  SHEAR  STRESS 
CALCULATIONS  FOR  EXPERIMENTAL  RUNS 


Table  B1 


Current  Velocity  Profiles  and  Shear  Stress 
Calculations  for  Experimental  Runs 


Set 

Run 

LP* 

■kie 

Y 

(cm) 

U1 

(cm/ sec) 

u*++ 

(cm/ sec ) 

(dynes/cm-' ) 

Buoy  C 

1 

b2 

57.9 

9.74 

45.7 

10.36 

30.5 

10.06 

15.2 

9.75 

6.1 

9.14 

0.24 

0.06 

Buoy  C 

2 

b2 

57.9 

17.43 

45.7 

18.20 

30.5 

15.18 

15.2 

14.33 

6.1 

14.23 

0.27 

0.08 

Buoy  C 

3 

c2 

57.9 

20.70 

45.7 

24.44 

30.5 

23.01 

15.2 

20.63 

6.1 

19.39 

0.83 

0.68 

Buoy  C 

4 

c2 

57.9 

31.91 

45.7 

34.41 

30.5 

32.10 

15.2 

28.96 

6.1 

25.27 

1.33 

1.75 

Lp  = Vertical  profile  location,  see  Figure  1 y. 


**  Y = Elevation  above  bed. 
t U = Horizontal  velocity. 

It  u*  = Friction  or  shear  velocity. 
I tq  = Bed  shear  stress. 


Fable  1U  (continued) 


Y 

U 

u* 

To 

Set 

Run 

LP 

(cm) 

(cm/sec) 

(cm/sec) 

(dynes/cm2 ) 

Buoy  C 

5 

E2 

57.9 

21.12 

45.7 

23.74 

30.5 

23.93 

15.2 

21.61 

6.1 

19.02 

1.17 

1.36 

Buoy  C 

6 

^2 

57.9 

9.72 

45.7 

10.29 

30.5 

10.10 

15.2 

9.70 

6.1 

9.08 

0.48 

0.23 

Block  15 

1 

57.9 

10.67 

45.7 

10.36 

30.5 

10.06 

15.2 

9.67 

6.1 

6.10 

0.14 

0.37 

Block  15 

2 

^2 

57.9 

15.82 

45.7 

16.58 

30.5 

16.58 

15.2 

7.1 

11.37 

0.78 

0.61 

Block  15 

3 

E2 

57.9 

21.12 

45.7 

23.74 

30.5 

23.93 

15.2 

21.61 

6.1 

19.57 

1.22 

1.49 

Block  15 

4 

E2 

57.9 

30.48 

45.7 

32.58 

30.5 

32.61 

15.2 

30.08 

6.1 

26.58 

1.59 

2.53 

Block  15 

5 

c? 

57.9 

9.75 

c 

45.7 

10.36 

30.5 

10.06 

15.2 

10.67 

6.1 

6.17 

0.91 

0.83 

B3 


Y U U*  T0 

Set  Run  lP  (cm)  (cm/sec)  (cm/sec)  (dynes/cm2) 


Block  27  1 Bo  54.9  13.11 

42.7  14.02 

30.5  13.41 

24.8  12.19 

21.3  11.58 

18.3  11.28 

15.2  11.57 

12.2  10.06 

9.1  9.75  0.24 

Block  27  2 B2  54.9  22.86 

42.7  23.77 

30.5  22.86 

24.8  22.56 

21.3  22.25 

18.3  21.95 

15.2  21.03 

12.2  20.12 

9.1  19.20  0.42 


Table  131  (continued) 


Y U u* 

Set  Run  Lp  (cm)  (cm/sec)  (cm/sec) 


Block  27  5 


Block  27  6 


Block  27  7 


Block  27  8 


Block  27  9 


24.8 

43.89 

21.3 

43.59 

18.3 

43.28 

15.2 

42.98 

12.2 

42.06 

9.1 

40.23 

54.9 

52.93 

42.7 

56.39 

30.5 

57.61 

24.8 

55.78 

21.3 

55.47 

18.3 

54.56 

15.2 

53.95 

12.2 

53.34 

9.1 

49.07 

54.9 

29.87 

52.7 

32.00 

30.5 

30.78 

24.8 

30.18 

21.3 

29.87 

18.3 

30.48 

15.2 

29.87 

12.2 

29.57 

9.1 

27.74 

54.9 

5.79 

42.7 

6.10 

30.5 

6.10 

24.8 

5.79 

21.3 

6.10 

18.3 

5.79 

15.2 

5.79 

12.2 

5.49 

9.1 

5.18 

54.9 

13.11 

42.7 

14.33 

30.5 

14.94 

24.8 

14.02 

21.3 

13.72 

18.3 

13.72 

1.44 


1.90 


0.84 


0.26 


B5 


(dynes/ cm2 ) 


2.07 


3.63 


0.71 


0.07 


Table  B1  (continued) 


Y 

U 

u* 

To 

Set 

Run 

LP 

(cm) 

(cm/sec) 

(cm/sec) 

(dynes/cm2 ) 

Block  D 

3 

b2 

54.9 

30.18 

42.7 

35.05 

30.5 

33.53 

24.8 

32.00 

21.3 

30.48 

18.3 

30.48 

15.2 

30. 18 

12.2 

29.87 

9.1 

29.57 

0.98 

0.96 

Buoy  D 

4 

b2 

54.9 

36.88 

42.7 

39.62 

30.5 

38.10 

24.8 

38.71 

21.3 

38.10 

18.3 

37.49 

15.2 

36.88 

12.2 

36.58 

9.1 

34.44 

1.09 

1.20 

Buoy  D 

5 

b2 

54.9 

42.37 

42.7 

43.89 

30.5 

44.50 

24.8 

43.59 

21.3 

42.98 

18.3 

42.37 

15.2 

42.06 

12.2 

41.45 

9.1 

37.19 

1.48 

2.20 

Buoy  D 

6 

b2 

54.9 

54.86 

42.7 

60.96 

30.5 

61.26 

24.8 

60.35 

41.3 

59.74 

18.3 

59.13 

15.2 

57.91 

12.2 

57.30 

9.1 

54.25 

3.0 

9.53 

B7 


A 


Y U U*  t0 

Set  Run  Lp  (cm)  (cm/sec)  (cm/sec)  (dynes/cm2) 


- 


, 


I 


Buoy  D 7 Bp  54.9 

42.7 

30.5 

24.8 

21.3 

18.3 

15.2 

12.2 

9.1 

Buoy  0 8 Bp  54.9 

42.7 

30.5 

24.8 

21.3 

18.3 

15.2 

12.2 

9.1 

Buoy  D 9 Bp  54.9 

42.7 

30.5 

24.8 

21.3 

18.3 

15.2 

12.2 

9.1 

Buoy  D 10  Bp  54.9 

42.7 

30.5 

24.8 

21.3 

18.3 

15.2 

12.2 

9.1 


20.12 

21.34 

20.73 

19.81 

19.51 

19.20 

18.29 

17.68 

17.07  1.19  1.42 

3.05 

3.35 

3.05 

3.05 

3.02 

3.05 

3.02 

2.99 

2.96  0.05  0.003 

14.63 

15.54 

15.85 
16.15 
15.54 

15.24 
14.94 
14.33 

13.41  0.80  0.63 

22.56 

23.47 

23.77 

23.77 

23.77 

23.47 

22.86 

22.25 

20.73  0.859  0.73 


3 


B8 


Table  B1  (concluded) 


Set  Run  Lp 


Y U u*  in 

(cm)  (cm/sec)  (cm/sec)  (dynes/cm2) 


Buoy  D 11  Bp  54.9 

42.7 
30.5 

24.8 

21.3 

18.3 

15.2 

12.2 

9.1 


30.48 

32.00 

32.00 

31.39 

31.09 

30.78 

30.48 

30.18 

29.87  0.875  0.76 


B9 


APPENDIX  C: 


EXPERIMENTAL  RUNS 


Figure  Cl.  Buoy  D dense  bed  erosional  runs  (visual  plot) 


23  cm /sec 


VaiN3DN03  WS1 


TIME  IN  HOURS  AFTER  VELOCITY  CHANGE 


s/un  2C 


Igure  C3.  Buoy  C dense  bed  erosional  runs 


Figure  C5.  Block  27  floe  bed  erosional  runs 


APPLNDI X D:  WASHLOAD  GRAIN-SIZE  ANALYSES 


Mashload  Grain-size  Analyses 
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Table  D1  (continued) 


Table  01  (concluded) 
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D5 


Equations  used  for  Reynolds  (Re)  and  Froude  (Fr)  Number  Calculations  were: 


Re  = 


ULp 

U 


Fr  « (gL) 


1/2 


where,  for  open  channel  flow  conditions: 

U = 

L = 61  cm  = average  water  depth 

p = 1.022  g/cm3  = density  of  salt  water  at  28  °/oo  and  25°C 

u = 9.5  x I0"3  g/sec  x cm  = dynamic  viscosity  of  salt  water  at 

28°/oo  and  25°C 

g = 980  cm/sec^  = acceleration  due  to  gravity 
and  for  offshore  bottom  flow  conditions: 

U = 75  cm/sec  = maximum  recorded  bottom  current  speed  (measured 
1.0  m above  bed) 

L = 300  cm  = approximate  height  of  bottom  currents  as  evidenced 
by  instantaneous  velocity  profiles  obtained  at  Buoy  C 
and  Buoy  D 
p * I .022  g/cm3 
li  s 9.5  x !0~3  g/sec  x cm 


g = 980  cm/sec 1 


E2 


I 


""" — 

Svmbo 1 

Mean  1 nq 

C' 

percentage  of  deposi table  sediment 

Co 

coefficient  of  cohesion 

cu 

remolded  shear  strength 

d 

particle  diameter 

i 

! 

orbital  diameter 

1 

e 

void  ratio 

Fr 

Froude  number 

9 

acceleration  due  to  gravity 

1 

] 

G 

specific  gravity 

i 

h 

water  depth 

H 

wave  height 

!P 

plasticity  Index 

jj 

k 

wave  number 

I 

I 

L 

reference  length 

ld 

deepwater  wavelength 

1 

Ls 

sample  location 

* 

LP 

profile  location 

Lw 

wave  length 

m 

slope 

1 

| 

Mz 

graphic  mean 

< 

Re 

Reynolds  number 

«h 

hydraulic  radius  or  wetted  perimeter 

S 

slope  of  the  energy  grade  line 

F2 

[ 

. jl 

Symbol 

SKj 

t 

T 

TSM 


Tt 

Ty 

U 

u 

^m 

U 

x 

u* 

Vs 

VT 

Vy 

w 


ws 

Wsw 


W 

Y 

Y 

Y 
Z 
Z 


T 

d 

w 

o 


Meaning 

inclusive  graphic  skewness 
time 

wave  period 

total  suspended  matter 

maximum  torque 

time  after  velocity  change 

horizontal  flow  speed 

mean  flow  spedd 

oscillatory  bottom  velocily 

horizontal  velocity  component 

friction  or  shear  velocity 

solid  particle  volume 

total  sediment  volume 

void  volume 

water  content 

I iquid  I imit 

plastic  I Imit 

weight  of  oven-dried  sediment 

weight  of  seawater 
fotal  weight  of  sediment  mass 
elevation  above  the  bed 
bulk  density 

specific  gravity  of  distilled  water  at  4°C 
water  depth 

roughness  element  = A = ks/x  In  Equation  13 


F3 


Symbol 

a 

H 

Y 

♦l 

P 

T 

O 

T 

ocr 

u 

v 

6 

0 


Mean  I nq 

constant 

constant 

un 1 1 weight  of  fluid 
Inclusive  graphic  standard  deviation 
water  density 
bed  shear 

critical  bed  shear 
dynamic  viscosity 
kinematic  viscosity 
laminar  sublayer 
rotation  angle 


F*1 


In  accordance  with  lot  tor  from  PAIN  RPC,  PAIN  AS1  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  m Library  of  Congress  MARC  format  is  reproduced 
be  1 ow . 
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